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ABSTRACT  '  _ _ _ 

j  Given  a  general  dynamic  equilibrium  formulation  of  a  time  staged  model,  JL  seek  conditions 
on  the  distribution  of  utility  functions  of  individuals  which  imply  the  model  is  equivalent  to  a 
mathematical  program.  In'-  ,  '1  ' 

Gorman  and  others  long  ago  have  observed  that  Engel  curves  of  average  consumption  as  a 
function  of  income  at  fixed  prices  are  remarkably  linear  over  a  broad  range  of  income  of  interest 

which  tapers  off  at  both  ends  of  this  range.  Wte  reproduce  this  phenomenon  by  assuming  - " 

that  a  general  polynomial  of  the  second  degree  has  enough  parameters  (coefficients)  to  globally 
represent  the  utility  functions  of  individual  consumers,  and  "(b)  the  distribution  of  utility  functions  _ 
that  individuals  have  is  independent  of  the  income  they  happen  to  have.  Wfachieve  the  latteFTry"  /  '  - 
assigning  values  to  the  parameters  of  the  utility  functions  by  a  random  drawing  with  replacement  f  — 
from  &  fepopulation  urn*’>containing  a  representative  sets  of  the  parameters.  We/then  derive  the 
functional  form  of  the  per  capita  demand  function  and  necessary  and  sufficient  conditions  for  its 
integrability.  >- 

Finally,  we  show;  in  the  context  of  the  time  staged  model,  that  when  the  population  is  not  too 
polarized  as  to  its  tastes  at  fixed  income  levels,  a  concave  objective  function  always  exists,  which 
maximized  subject  to  the  physical  flow  constraints,  implies  the  equilibrium  conditions.  • 
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Introduction 


Given  a  general  dynamic  equilibrium  model,  a  long  standing  problem  is  finding  conditions  on 
the  distribution  of  utility  functions  of  individuals  which  guarantee  that  the  model  is  equivalent  to 
the  problem  of  maximizing  a  concave  objective  function  of  the  aggregate  consumption  variables 
subject  to  the  physical  constraints  of  the  system.  When  the  two  are  equivalent,  the  powerful 
software  of  mathematical  programming  can  be  applied  to  efficiently  solve  large  scale  equilibrium 
problems. 

Our  approach  differs  from  past  ones  by  the  way  we  assign  values  to  the  set  of  a  parameters  of 
the  utility  function  of  individuals.  We  first  assume  the  distribution  is  independent  of  the  income 
individuals  happen  to  have.  We  achieve  this  independence  by  placing  representative  vectors  of 
parameters  values  in  an  “urn”  and  assigning  them  to  each  individual  at  an  income  level  by  a 
random  drawing  with  replacement  from  the  urn. 

1  The  authors  wish  to  thank  Kenneth  Arrow,  Gerard  Debreu,  Robert  Dorfman,  Dale  Jorgenson 
and  Lawrence  J.  Lau  for  their  helpful  comments. 
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In  the  context  of  a  time-staged  equilibrium  model,  we  are  therefore  not  interested  in  deriving 
the  functional  form  of  the  per  capita  demand  function  of  some  future  time  period,  since  it  is  now 
a  random  function,  but  we  are  interested  in  expected  demand  as  a  function  of  prices  and  per 
capita  income.  By  the  law  of  large  numbers,  the  expected  per  capita  demand  will  differ  from  actual 
future  per  capita  demand  insignificantly  for  a  population  the  size  of  the  U.S.  Accordingly  we  will 
be  seeking  conditions  on  the  distribution  of  representative  utility  function  parameters  in  the  urn 
which  imply  that  these  expected  per  capita  demand  functions  are  integrable. 

An  early  controversy  which  arose  in  theoretical  work  on  consumer  demand  concerned  the 
shape  of  Engel  curves  which  express  at  fixed  prices  average  consumption  of  a  particular  category  of 
goods  (such  as  food)  as  a  function  of  income  level.  Survey  data  strongly  suggest  that  these  curves 
are  nearly  linear  over  a  broad  range  of  income  levels,  see  Graphs  1-8  at  the  end  of  this  paper. 
Some  authors  have  conjectured  that  his  must  be  true  because  the  underlying  demand  functions  of 
individuals  must  be  linear,  or  nearly  so,  at  fixed  prices  over  a  broad  range  of  individual  income 
available  for  consumption. 

For  example,  Gorman  in  his  1953  paper  [9]  remarks:  “A  great  deal  of  work  has  been  done  on 
Engel  curves  particularly  by  Allen  and  Bowley  and  Houthakker.  The  work  of  Allen  and  Bowley 
was  based  on  the  assumption  that  the  classical  Engel  curves  for  different  individuals  at  the  same 
prices  were  parallel  straight  lines,  but  this  has  been  rejected  in  work  of  Houthakker  in  favor  of  a 
doubly  logarithmic  form.  However,  the  earlier  assumption  fits  the  data  remarkably  well.”  [22,23], 

Gorman,  in  the  article  cited  above,  finally  showed  conclusively  that  average  linear  demand 
could  only  happen  if  the  underlying  individual  demand  functions  are  linear  in  income  at  fixed 
prices. 

Since  linearity  of  individual  demand  functions  in  income  at  fixed  prices  was  considered  unlikely 
for  theoretical  reasons,  researchers  have  concentrated  on  using  the  logarithmic  functional  form 
following  the  traditional  approach  to  demand  analysis  found  in  the  pioneering  work  of  Henry 
Shuitz,  Richard  Stone,  and  Herman  Wold,  [24,25,26]. 

The  logarithmic  utility  functions  were  generalized  by  Christianson,  Jorgenson,  and  Lau  [28,30] 
who  proposed  that  utility  be  approximated  by  a  function  which  is  quadratic  in  the  logs  of  the 
sumption  variables.  Certain  difficulties  in  translating  from  individual  demand  to  aggregate 
demand  which  were  present  in  their  early  work  were  overcome  by  a  more  general  theory  which 
the  authors  Jorgenson  Lau  and  Stoker  call  “exact  aggregation”,  see  [11]  They  consider  a  model 
in  which  each  individual  has  a  trans-log  utility  function  which  depends  on  “attributes”  of  the 
individual  as  well  as  consumption  itself. 

By  placing  restrictions  on  the  way  that  these  attributes  enter  into  the  utility  function,  they 
were  able  to  find  a  model  in  which  demand  exactly  aggregates  thereby  deriving  per  capita  and 
aggregate  demand  as  a  function  of  prices  and  a  certain  class  of  symmetric  statistics  of  the  income 
distribution,  which  includes  more  than  just  the  mean  or  total  aggregate  income.  In  this  way  they 
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are  able  to  obtain  an  exact  aggregation  result  without  appealing  to  linearity  in  individual  or  average 
demand  at  various  income  levels. 

Our  approach  is  motivated  theoretically  by  the  observation  that  when  the  utility  functions 
of  individuals  are  polynomial  expressions  of  the  second  degree  in  the  consumption  variables,  the 
demand  function  of  each  individual  will  be  linear  over  a  broad  range  of  income  at  fixed  prices 
tapering  off  to  zero  at  the  poverty  end  as  more  and  more  individuals,  in  order  to  maximize  their 
utility,  have  to  set  more  and  more  components  of  their  consumption  vector  to  zero,  and  leveling 
off  at  the  high  income  end  as  more  and  more  individuals  have  income  sufficient  to  purchase  their 
“satiation”  vector.  Consistent  with  empirical  observation,  this  implies  that  expected  demand 
functions  (Engel  curves)  are  linear  over  a  broad  range  of  income  (used  for  consumption),  and  this 
in  turn  implies  that  the  expected  per  capita  demand  is  linear  in  per  capita  income  over  a  broad 
range  of  per  capita  income. 

There  is  no  obvious  reason  why  these  steps  should  lead  to  an  expected  per  capita  demand 
function  for  a  given  time  period  t  that  is  integrable;  or  if  integrable  within  a  period  why  this  should 
imply,  in  the  context  of  a  time-staged  equilibrium  model,  that  one  can  replace  the  dual  price 
constraints  by  an  objective  function  which  can  be  maximized  subject  to  just  the  primal  system 
of  physical-flow  constraints  to  obtain  the  equilibrium  solution.  Indeed  from  Arrow’s  Impossibility 
Theorem,  we  know  that  seeking  an  aggregate  utility  function  for  the  economy  over  time  could  be 
a  futile  quest  since  in  general  it  need  not  exist  [l,  5,  18,  27). 

Conditions  for  integrability  of  a  demand  function  have  been  given  by  Slutsky,  see  Varian  [19], 
We  derive  a  necessary  and  sufficient  condition  for  integrability  based  on  the  derived  form  of  the 
expected  aggregate  demand  function  for  period  t. 

It  all  depends  on  how  the  parameters  of  the  utility  functions  of  individuals  in  urn  are  dis¬ 
tributed.  One  measure  of  how  much  utility  functions  £7’  differ  from  one  another  is  to  compare  H', 
where  Hl  is  the  price  cross-effect  matrix  of  the  «-th  utility  function  in  the  urn  with  H  any  positive 
definite  matrix  used  for  comparison  (such  as  average  H')  by  forming  H'  =  tf-1/2#*#-1/2  .  We 
prove,  in  a  worse  case  scenario,  that  if  p,  <  3  +  Vs  =  5.83  for  every  where  pt  ,  is  the  ratio  of  the 
largest  to  smallest  eigenvalue  of  H' ,  then  the  expected  per  capita  demand  function  is  integrable. 
Moreover,  the  greater  the  variability  of  orientations  of  the  axes  of  the  ellipsoids  p‘ H'p  =  constant 
the  higher  is  the  bound  for  pi. 

Outline  of  the  Paper 

in  Part  I,  we  derive  the  form  of  the  expected  per-capita  demand  function  for  each  discrete  time 
period  t;  in  Parts  II  and  III  we  derive  necessary  and  sufficient  conditions  for  its  integrability;  in  Part 
IV,  we  present  the  time-staged  dynamic  equilibrium  model  for  consumers  and  producers /investors. 
In  Part  V,  we  develop  a  first-order  approximation  to  the  expected  aggregate  demand  and  utility 
functions  for  a  period  which  allows  one  to  explicitly  express  these  functions  and  estimate  their 
parameters.  Finally  we  estimate  the  parameters  of  the  expected  per  capita  demand  function  and 


the  per  capita  utility  function  using  survey  data  and  test  the  theory  empirically  by  using  it  to 
predict  time-series  data  of  per  capita  consumption  of  various  iterrk.  when  prices  and  per  capita 
income  are  given. 

Even  when  the  latter  exists,  if  the  numeraire  for  normalizing  period  t  prices  is  not  suitably 
chosen,  we  prove  in  the  context  of  time-staged  equilibrium  model  that  a  utility  function  that  drives 
the  economy  over  time  need  not  exist.  However,  when  a  numeraire  is  suitably  chosen,  the  dynamic 
equilibrium  model  is  integrable. 

Prices  normalized  by  this  numeraire  differ  very  little  in  practice  from  those  obtained  by  scaling 
them  so  that  their  average  price  is  unity,  see  last  column  of  Table  4.  Thus,  from  the  viewpoint 
of  the  investor,  these  prices  are  equally  acceptable  for  calculating  the  rates  of  return  of  various 
investment  possibilities.  If  it  is,  then  it  is  possible  to  restate  the  dynamic  equilibrium  problem  as 
a  mathematical  program  and  to  use  non-linear  programming  software  like  MINOS  to  optimize  the 
primal  system  [15]. 

To  test  the  theory,  fits  were  made  to  survey  data  (see  Table  I)  and  tested  by  predicting  the 
consumption  pattern  of  final  consumers  for  the  years  1961  to  1982  as  a  function  of  prices  and  per 
capita  income  used  for  consumption  (see  Graphs  9  to  16).  We  also  report  on  an  experiment  that 
suggests  the  approximation  may  be  a  very  good  one  even  if  the  set  of  utility  functions  U *  in  the 
urn  differ  markedly  from  one  another. 

Motivation  of  our  Research 

Our  presentation  here  arose  out  of  our  efforts  beginning  in  1975  to  build  a  macro-economic 
model  of  U.S.  to  assess  the  long-  term  effects  of  modernization,  innovation,  foreign  competition, 
energy  prices,  and  conservation  on  the  growth  of  various  economic  sectors,  GNP,  and  per  capita 
income.  PILOT,  is  a  multi-time  period  model,  quite  large,  with  a  data  base  of  over  70,000  tech¬ 
nological  coefficients,  [4].  Its  principal  weakness,  as  we  see  it,  lies  not  with  the  numerical  results 
from  various  scenario  runs  (the  physical  growth  of  the  economy  appears  to  be  quite  reasonable), 
but  with  our  inability  to  justify  the  aggregate  utility  function  which  we  had  devised  to  increase  the 
standard  of  living  and  had  been  using  as  the  driver  as  though  the  U.S.  were  a  planned  economy. 
To  be  more  precise,  the  partial  derivatives  of  this  somewhat  arbitrarily  chosen  utility  function, 
interpreted  as  prices,  implied  a  behavioral  response  of  final  consumers  to  prices  and  a  behavioral 
response  of  produce rs/investors  also  to  prices,  which  are  almost  certain  to  be  out  of  kilter  with 
what  their  observed  behaviors  would  be  in  the  real  world. 

These  considerations  led  us  recently  to  reverse  the  process  and  to  reformulate  PILOT  along 
classical  economic  lines  as  a  dynamic  equilibrium  model  that  satisfies  the  behavioral  responses 
of  final  consumers  to  prices  given  their  income  for  consumption,  and  the  responses  of  produc¬ 
ers/investors  in  choosing  activities  that  yield  at  least  a  minimum  rate  of  return.  From  a  mathe¬ 
matical  point  of  view,  this  is  not  a  dramatic  change,  since  the  set  of  Kuhn- Tucker  conditions  [12] 
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that  must  be  satisfied  at  the  point  which  maximizes  >  he  utility  function,  are  quite  analogous  to  the 
Arrow-Debreu  conditions  [1,5]  that  must  be  satisfied  at  an  equilibrium  point.  The  main  difference 
is  that  the  general  equilibrium  problem  belongs  mathematically  to  the  more  general  class  of  com¬ 
plementarity  problems  which  require  for  solution  combining  the  dual  system  of  price  constraints 
along  with  the  primal  system  of  physical  constraints  in  one  big  simultaneous  system,  [3,  6,  8,  13, 
17].  For  a  problem  the  size  of  our  PILOT  model,  this  combined  system  is  too  large  to  solve  directly 
and  much  of  our  research  has  been  concerned  with  finding  efficient  ways  to  use  mathematical  pro¬ 
gramming  software,  like  MINOS  [15],  to  solve  the  system,  including  ways  to  decompose  the  system 
into  smaller  problems  and  to  use  their  solutions  to  iteratively  converge  to  an  equilibrium. 

A  fundamental  question  that  has  concerned  us  is  reconciling  the  “prescriptive”  (normative) 
view  of  the  initial  formulation  with  the  “descriptive”  (behavioral)  view  of  the  dynamic  equilibrium 
formulation.  To  be  precise,  do  the  realistic  behavioral  assumptions  of  an  equilibrium  model  serve  as 
a  driver  promoting  reasonable  growth  and  well  being  of  the  economy  when  it  has  the  potential  for 
such  growth?  We  will  see  later  that  the  dynamic  equilibrium  formulation  implies  under  reasonable 
conditions  an  objective  for  the  economy  the  form  of  which  makes  it  evident  why  growth  in  the 
economy  will  occur  if  it  has  potential  for  growth. 
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PART  I: 

DERIVING  THE  PER  CAPITA  DEMAND  FUNCTION  FOR  PERIOD  t 


Notation 

For  discussion  involving  a  fixed  period,  we  will  usually  omit  the  time  subscript.  Thus  the  con¬ 
sumption  vector  of  the  j- th  consumer  in  period  t  is  denoted  by  X3  =  X} ;  his  budget  or  personal 
income  for  consumption  is  denoted  by  I  =  It  measured  in  period  t  undiscounted  dollars.  Expected 
per  capita  consumption,  income,  and  utility  we  denoted  by  Xt,It,Ut  and  the  corresponding  ag¬ 
gregates  by  bold  face  X*  =  Pt  •  Xt,  It  =  Pf  It,  U t  =  Pt  -Ut  where  Pt  is  the  size  of  the  population 
in  period  t. 

We  will  use  the  symbol  I  to  denote  the  identity  matrix  to  avoid  confusion  with  I  and  I  which 
refer  to  income.  The  inner  product  of  a  column  vector  v  with  itself  will  be  denoted  by  v2  =  v'v.  The 
symbols  a, 0,6,  A  denote  scalar  constants.  L.H.S  and  R.H.S.  are  abbreviations  for  left  hand  side 
and  right  hand  side  of  an  equation  or  an  inequality  relation.  The  symbol  =  means  approximately 
equal. 

Utility  Function  of  Individuals 

The  first  assumption  we  make  is  that  each  final  consumer  has  a  utility  function  that  is  quadratic 
in  the  consumption  variables  and  that  he  chooses  his  consumption  vector  by  maximizing  his  utility 
function  subject  to  his  budget  constraint.  This  specific  functional  form  may  be  viewed  as  a  second- 
order  global  approximation  to  whatever  may  be  his  true  utility  function. 

Assumption  1.  Individual  j  in  period  t  has  a  utility  function  U3(X),  measuring  the  value  j 
attaches  to  having  a  consumption  vector  X,  which  can  be  represented  by  a  general  quadratic 
function  of  the  form 

U3(X)  =  2 {M3S3)'X  -  [X)'M3{X)  +  Constant ,  (1.1.0) 

where  vector  S3  >  0  and  matrix  M3  is  symmetric  and  positive  definite,  hence  non-singular. 

Without  loss  of  generality,  we  may  rescale  the  matrices  M3  so  that,  letting  e'  =  (1,1,...,  1), 
their  inverses  H3  have  the  property 

=  =  1,  H^iM3)-1.  (1.1.1) 

k  l 

Letting  Constant y  =  ~(S3)'M3 S3 , 

U3{X)  =  -( 5 3  -  X)'M3(S3  -X)<0,  for  all  S3  -  X  ?  0.  (1.1.2) 

It  is  easy  to  see  U 3  ( X )  in  (1.1.2)  is  unconditionally  maximized  when  X  —  S3 .  Therefore,  it  is  natural 
to  assume  that  S  is  strictly  positive  and  to  refer  to  S3  as  the  “satiation”  vector  of  individual  j.  We 


view  the  income  /  of  an  individual  as  an  “authorization”  to  expend  up  to  that  amount  for  actual 
consumption.  Should  it  happen  that  an  individual’s  budget  I  >  p'S3 ,  he  maximizes  his  utility  by 
buying  his  satiation  vector  S3 .  The  unexpended  amount  I  -  p'S3 ,  in  this  case,  is  not  used. 

Form  of  the  individual  demand  function.  The  budget  constraint  and  nonnegativity  constraint 
for  individual  j  in  period  t  at  fixed  prices  p  are: 

p'X<I ,  X>0 ,  p>  0.  (1.2.0) 

Subject  to  (1.2.0),  j  maximizes  his  utility  function  U3(X).  The  vector  X  which  maximizes  U3  (X) 
subject  to  a  budget  constraint  is  denoted  by  X3 .  This  is  a  special  case  of  a  quadratic  programming 
problem  [3],  We  distinguish  three  cases:  The  high  income  case,  the  standard  case  that  defines 
the  range  of  income  of  interest,  and  the  low  income  case.  In  the  high  income  case,  I  >  p'S3  and 
individual  j  maximizes  his  utility  by  buying  his  satiation  vector,  i.e.,  X3  =  S3 .  Otherwise,  the 
budget  constraint  is  tight  and  the  procedure  begins  by  forming  a  Lagrangian  and  setting  its  partials 
to  zero, 

d[U3{X)  -  2\(p'X)\/dX  =  0,  (1.3.0) 

which  is  then  solved  to  determine  X3  as  a  function  of  A;  the  expression  for  X3  as  a  vector  function 
of  A  is  then  substituted  into  (1.2.0)  with  the  budget  tight  and  solved  for  A.  If  X3  >  0,  then  this 
is  the  standard  case.  The  low  case  occurs,  by  definition,  when  the  budget  1  <  /*  is  so  low  that  j 
maximizes  his  utility  on  the  boundary  of  the  non-negative  orthant  X  >  0  by  setting  to  zero  one  or 
more  components  of  X. 

We  will  loosely  refer  to  the  income  levels  between  /*  =  max  /*  and  I**  —  min  p'S3  as  the 
range  of  income  of  interest.  It  is  the  range  of  income  for  consumption  in  which  no  individual  j 
maximizes  his  utility  on  the  boundary  of  the  orthant  X  >  0  or  has  sufficient  income  to  buy  his 
satiation  vector  S3 .  This  range  depends  on  the  prices  p.  We  think  of  this  range  as  very  broad,  I* 
representing  extreme  poverty  and  /**  as  being  very  rich.  In  the  context  of  the  full  model  (Part 
IV),  it  is  possible  to  have  all  individuals  receive  at  least  a  minimum  fixed  consumption  vector,  so 
that  X3  (k)  =  0  should  not  be  interpreted  as  j  going  without  food,  for  example.  Income,  in  this 
context,  means  income  for  purchases  above  this  floor. 

For  the  “standard”  case,  we  substitute  the  quadratic  expression  for  U3{X)  given  by  (1.1.0) 
into  (1.3.0)  and  differentiate  partially: 

M3{S3  -  X3)  =  A-p.  (1.4.0) 

Solving  for  S3  -  X3 : 

S3  -  X3  =  A  •  H3p,  where  H3  =  {M3)~\  (1.5.0) 

Note  that  the  inverses  H3  =  (M;)-1  exist  and  are  also  symmetric  and  positive  definite.  We  can 
now  use  the  tight  budget  constraint  to  determine  A.  Multiplying  (1.5.0)  by  p'  on  the  left  and 


I 


I 

| 

& 

$ 

1 

l 


a 


I 


setting  =  I,  we  can  solve  for  A  and  substitute  back  into  (1.5.0).  This  yields  (1.6.0)  below,  the 
demand  function  of  individual  j  as  a  function  of  prices  and  income  I  where  A,  equal  to  the  term  in 
parentheses,  is  positive  if  the  budget  I  is  less  than  the  income  p'S 3  required  to  buy  the  satiation 
vector  S3 .  Note  p'H3p  >  0  follows  from  H3  being  positive  definite.  It  is  easy  to  see  that  (1.6.0) 
satisfies  two  properties  of  demand  functions  which  hold  when  the  budget  constraint  is  tight:  (a) 
p' X 3  —  I  and  (b),  X1  remains  invariant  if  we  rescale  prices  p  and  income  7  proportionally. 

Theorem  1.1.  The  demand  function  of  individual  j,  for  fixed  prices  p,  is  a  piecewise  linear  function 
of  consumption  income  I.  (i)  For  a  range  that  includes  the  “range  of  income  of  interest”,  the  demand 
function  is  given  by 


si  -  *  = 


/;  <  /  <  /;*  =  Ps 3 


(1.6.0) 


(ii)  For  I  <  /*,  one  or  more  components  X3(k)  =0.  As  I  — ►  0,  X3  — *  0.  (Hi)  For  I  >  7*‘  =  p'S3 , 
the  budget  constraint  is  slack  and  X3  =  S3  for  all  I  >  /** . 


Proof.  The  theorem  is  a  restatement  of  a  well  known  result  from  quadratic  programming.  We 
therefore  omit  details  and  sketch  only  important  steps  in  the  proof.  The  Kuhn-Tucker  conditions 
for  testing  whether  X  =  X(I)  is  optimum  consists  in  partitioning  the  inequalities  generated  by 
taking  partial  derivatives  into  a  set  fl  that  are  tight  under  the  current  solution  and  Cl,  those 
that  are  not.  These  are  in  1  —  1  correspondence  with  components  Xk  >  0  and  =  0.  Call 
the  corresponding  index  sets  k  also  Cl  and  Cl.  Suppose  for  two  income  levels  7  =  a  and  I  =  (3 
the  corresponding  index  sets  fl  are  the  same  for  the  two  optimum  solutions.  Then  the  convex 
combination  X  =  AX(a)  +  pX(fl),  I  =  Aa  +  p/?,  is  optimal  for  all  A  +  p  =  1,  (A,p)  >  0.  It  follows 
that  X  is  linear  in  I  over  segments  where  fl  sets  are  the  same,  there  being  at  most  one  segment  for 
each  such  fl.  For  other  fl  there  may  be  only  one  I.  The  former  are  the  broken  line  segments  and 
the  latter  are  the  values  of  I  corresponding  to  the  break  points.  Since  the  number  of  different  fl  is 
finite,  there  are  a  finite  number  of  breakpoints  and  broken  line  segments.  | 

Expected  Demand  Function  for  Individuals  at  the  Same  Income  Level.  Our  next  step  is 
to  derive  the  functional  form  of  the  Engel  curves,  namely  the  average  demand  for  some  item  k  for 
all  individuals  with  the  same  income  level  7  as  a  function  of  7  and  prices  p.  For  this  purpose,  as 
we  already  noted,  we  assume  that  the  utility  functions  which  people  have  are  independent  of  the 
particular  income  which  they  happen  to  have.  Suppose  there  are  n/  individuals  j  at  income  level  7. 
We  need  to  have  a  way  to  assign  the  parameters  (S3 ,  H3)  to  their  utility  functions  independent  of 
their  income  7,  next  find  X3  =  X  which  maximizes  their  utility  subject  to  their  budget  constraint 
7,  and  then  average  X3  to  compute  the  average  demand  vector  as  a  function  of  7  and  fixed  prices 


Population  Urn.  We  can  achieve  this  independence  of  income  by  selecting  in  the  parameter 
space  ( S,H )  a  representative  set  of  n  possible  points  (S',  H')  and  writing  them  as  labels  on  n 
balls,  t  =  (1, . . .  ,n),  and  placing  the  balls  in  an  imaginary  urn  which  we  call  the  population  urn.  In 
making  a  random  drawing  with  replacement  from  the  urn,  we  assume  each  (S',  H')  is  equally  likely 
to  be  drawn  and  assigned  to  individual  j.  If  it  is  desired  to  make  it  more  likely  to  choose  certain 
(S',H'),  this  could  be  done  by  weighing  the  distribution  of  (S',  H')  in  the  urn  or  by  replication 
of  certain  of  the  balls.  For  convenience  we  have  made  this  distribution  discrete  and  n  finite,  but  a 
general  probability  measure  could  be  used  instead. 


Assumption  2,  the  Population  Urn  Assumption.  The  utility  function  which  individual 
j  happens  to  have,  is  uncorrelated  with  I,  his  personal  income  for  consumption.  We  achieve 
this  by  assigning  to  j  a  utility  function  with  parameters  (S},  H°)  =  (S',  H')  where  (S' ,H')  is 
randomly  drawn  with  replacement  from  an  “urn”;  we  also  assume  the  distribution  in  the  urn  of 
the  “ satiation-level ”  parameters  S'  is  independent  of  the  distribution  of  the  price  “ cross-  effect” 
parameters  H'  =  (M*)-1. 

For  any  of  the  n/  individuals  j,  with  income  I,  we  denote  their  expected  S}  and  X 3  by 
£  S*  =  S1  and  £  X3  =  X1  respectively.  Our  assumption  that  the  distribution  of  utility  functions 
does  not  depend  on  /  implies  that  this  expected  S1  is  the  same  for  all  I  and  therefore  S1  -  S 
where  S  is  the  arithmetic  mean  of  the  S'  in  the  population  urn. 

A  special  symbol  is  used  to  denote  the  arithmetic  mean  of  S'  in  the  population  urn.  It 
has  the  same  value  as  the  expected  value  of  a  single  random  drawing  of  S'  from  the  population 
urn.  Likewise,  we  denote  the  arithmetic  mean  of  any  function  rp(S',H')  for  »  =  l,...,n  in  the 
population  urn  by  £itp(S' ,  H')  =  (1/n)  £3”  rp(S',Hx).  For  the  nj  individuals  j  that  have  an  income 
/  in  the  range  of  income  of  interest,  max  /*  =  /*</</**  =  min  /**,  we  know  that  (1.6.0)  holds 
no  matter  what  (S',  IP)  has  been  assigned  to  j  from  the  urn.  Therefore  taking  the  expectation  of 
(1.6.0)  we  obtain  (1.7.1)  below. 

Theorem  1.2.  For  fixed  prices  p  =  pt  and  for  all  individuals  j  whose  income  level  I  satisfies 
maxi ’*  =  /*  <  I  <  /**  =  min p'S*  for  (S',H')  in  the  population  urn,  their  expected  consumption 
vector  X1  is  a  linear  function  of  I: 


=  (pS-I)£i\-~Hi)p,  S  =  £iSi  . 


(1.7.1) 

(1.7.2) 


Proof.  Our  assumption  that  the  parameters  ( S3,H 3)  assigned  for  individual  j  is  chosen  by  a 
random  drawing  from  the  population  urn  independent  of  income  I  means,  for  fixed  prices,  that  the 
expected  values  of  the  terms  of  (1.6.0)  are  those  of  (1.7.1);  in  particular  the  factor  multiplying  I 
in  the  second  term  does  not  depend  on  I.  Therefore  (1,7.1),  at  fixed  prices,  states  that  the  Engel 


curves  which  express  each  component  k  of  X1  as  a  function  of  I  are  linear  in  /  at  fixed  prices  in 
the  range  /*  <  /  <  /**.  Our  additional  assumption  about  the  independence  of  the  distributions  of 
S'  and  H'  in  the  urn  implies  (1.7.2)  because 


e'{p's')[yw~PH']  p  =  £t(p,s'){£'7W~PH']p = {p^){£' 


p'H'p 


H']p 


(1.7.3) 


We  denote  the  difference  between  the  expected  demand  at  incomes  pS  and  pS  —  1  by 

G^)  =  e\~-Hi\p  (1.8.0) 

It  follows  that  the  difference  between  the  expected  demand  at  saturation  incomes  p'S  and  I  is 
J  ■  G(p),  where  J  =  pS  -  I  is  the  additional  income  required  to  reach  an  income  sufficient  to  buy 
the  satiation  vector.  Substituting  (1.8.0)  into  (1.7.2),  a  more  complete  statement  of  Theorem  1.2 
is 


Theorem  1.3.  For  fixed  prices  p  =  pt,  then  (i),  for  all  individuals  j  whose  income  level  I  is  in  the 
range  of  income  of  interest,  their  expected  consumption  vector  X1  is  a  linear  function  of  I: 

S  -  X  =  ( p'S  -  I)  ■  G{p),  where  max  /*  =  /*</</**=  minp'5*  .  (1.9.0) 

(ii)  .is  I  decreases  below  /*,  it  is  more  and  more  likely  individuals  j  with  this  income  will  maximize 
their  utility  by  setting  components  of  XJ  to  zero  and  £ X3  =  X1  —*  0.  (Hi)  As  I  increases  beyond 
/**  =  min  p'S' ,  it  is  more  and  more  likely  individuals  j  will  maximize  their  utility  by  buying  their 
satiation  vectors  and  EX3  =  X1  will  level  off  to  S. 


Com  ent.  In  Part  V,  we  will  present  the  empirical  data  of  average  consumption  as  a  function 
of  income  level  for  eight  broad  consumer  categories  Food,  Clothing,  Housing,  Household  Operation, 
Transportation,  Recreation,  Personal  Care,  and  All  Other.  We  will  see  that  the  form  of  these 
Engel  Curves  is  generally  consistent  with  this  theoretical  result  that  the  expected  demand  vector 
of  individuals  with  same  consumption  income  /  is  a  linear  function  of  I  over  the  “range  of  income 
of  practical  interest,”  namely  maxi *  <  J  <  min  p'S'. 


Per  capita  demand  as  a  function  of  per  capita  income.  Assuming  fixed  prices  in  period  t, 
the  expected  per  capita  deman  J  function  is  derived  from  the  expected  demand  of  persons  at  various 
income  levels  by  a  convolution  with  the  income  distribution.  Let  X1  (k)  =  C*(/)  be  the  expected 
personal  consumption  per  year  of  item  k  by  individuals  at  the  same  income  level  I.  Let  </>(/)  be 
the  income  distribution.  The  per  capita  consumption  of  item  k  and  corresponding  expected  per 
capita  income  per  year  are  given  by 


x{k)=  r  <t>{i)-ck{i)di, 

Ji- o 


(1.10.1) 

(1.10.2) 
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where  X(k)  denotes  component  k  of  X.  The  symbol  for  expected  per  capita  consumption  vector 
X  is  to  be  distinguished  from  X1  which  is  the  expected  consumption  vector  of  individuals  j  whose 
income  level  is  7.  It  would  appear  that  the  correspondence  between  X  and  7  depends  on  the 
distribution  of  <f>\  in  fact  it  does  not  for  a  broad  range  of  7. 

Theorem  1.4.  If  Ck(I)  is  a  linear  function  of  I,  then  independent  of  the  distribution  of  income 
<f>(I),  X(k)  =  C'fc(T),  i.e.,  expected  per  capita  consumption  of  item  k  is  a  linear  function  of  per 
capita  income  I  and  this  linear  function  is  Ck(T). 

Proof.  Let  C*(7)  =  o  +  bl.  Substituting  into  (1.10.1)  and  noting  /  <t>(l)dl  =  1,  yields  X(k)  — 
a  +  bI=Ck(T).  | 

Comment.  The  hypothesis  that  Ck(I)  is  a  linear  function  of  I  is  only  true  by  Theorems  1.2  and 
1.3,  for  a  restricted  range  of  income  which  we  have  referred  to  as  the  “broad”  range  of  income  of 
interest,  max  If  =  I*  <  I  <  /**  =  min  pS*  for  (S‘ ,  77')  in  the  population  um.  Let  us  assume 
the  distribution  of  income  for  consumption  in  the  population  is  above  the  extreme  poverty  level 
7*  =  max  If  and  below  being  very  rich,  /**  =  min  pS' .  As  time  goes  by,  per  capita  income  7  will 
change  (and  likely  increase)  and  the  distribution  of  income  <j>(I)  about  7  will  change.  As  long  as 
people  at  the  same  income  level  at  the  same  prices  buy  in  the  same  way  in  the  future  and  the 
income  distribution  change  is  not  so  drastic  that  some  j  have  income  7  below  their  7*  or  above 
their  pS3 ,  Theorem  1.4  states  that  X(k)  is  the  same  linear  function  of  7  as  X1  (k)  is  of  7  for  some 
range  7*  <  7  <  7**.  Therefore,  we  have  established: 

Theorem  1.5.  For  Sxed  prices  p  =  pt,  (i)  expected  per  capita  consumption  X  is  a  linear  function 
of  per  capita  income  I  for  a  certain  range  of  income  I;  namely: 

S  -  X  =  (p'S  -  7)  •  G(p)  ,  r  <  7  <  7**  ,  (1.11.1) 

where  G(p)  =  H(p)  •  p  and  II (p)  —  £i(p'/7’p)_1  •  77*  is  a  symmetric  positive-definite  matrix  whose 
elements  depend  on  p  and  not  on  I. 

(ii)  If  I  <  7*,  then  it  is  more  and  more  likely  that  some  individual’s  income  I  <  I*,  and  these 
j  will  maximize  their  utility  by  setting  some  components  of  X3  to  zero;  as  7  — »  0,  X  — *  0.  (Hi) 
If  7  >  /**,  then  it  is  more  and  more  likely  that  some  individual’s  income  I  >  p'S3;  these  j  will 
maximize  their  utility  by  buying  their  satiation  vectors  S3;  for  sufficiently  high  7,  X  will  level  off 
to  S. 

Note  especially  that  our  predicted  result  of  linearity  does  not  depend  on  the  shape  of  the 
income  distribution  or  how  this  changes  as  per  capita  income  I  increases  in  the  future  providing  it 
does  not  rise  so  high  that  some  j  have  income  7  such  tht  p'S3  <  7  or  decreases  so  low  that  some 
If  >  I.  Implicit,  of  course,  is  the  assumption  that  tastes  don’t  change  with  time.  If  they  do,  there 
is  no  problem  adjusting  the  model  for  trends  in  taste. 
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PART  II: 

NECESSARY  AND  SUFFICIENT  CONDITION  FOR  INTEGRALITY 
OF  THE  EXPECTED  PER-CAPITA  DEMAND  FUNCTION 


Our  immediate  objective  is  to  derive  a  necessary  and  sufficient  condition  that  an  inverse  ex¬ 
pected  per  capita  demand  function  and  utility  function  exist  for  period  t  when  the  expected  demand 
function  given  by  (1.11.1).  The  condition  for  existence  is  closely  related  to  one  given  by  Slutsky, 
see  Varian  [19].  Our  proof  is  specific  to  the  demand  function  given  by  (1.11.1).  It  is  therefore  as¬ 
sumed  in  the  discussion  that  follows  that  per  capita  income  is  in  the  “range  of  income  of  interest” 

r  <I<  7**. 

Definition  of  a  Utility  Function.  In  order  for  U (X)  to  qualify  as  a  per  capita  utility  function, 
we  let  v  =  S  —  X  and  require  Z(v)  =  ~U(X)  to  be  a  convex  function  twice  differentiable  (except 
possibly  at  v  =  0)  which  attains  a  minimum  subject  to  a  budget  constraint  p'X  =  7,  or  equivalently 
p'v  =  J ,  where  J  =  p'S  —  7,  at  a  unique  finite  point  v  =  v*.  Therefore  v*  is  a  function  of  p  and  J 
which  we  denote  by  S  —  X  =  v(p,  J)  and  call  the  latter  the  demand  function  associated  with  the 
utility  function  U(X). 

Conversely,  if  we  are  given  a  demand  function  S  -  X  =  v(p,  J),  such  as 

S-X  =  v  =  J  Sip'Wp)-1  •  Wp  where  J  =  p'S-7,  (2.1.0) 

t 

we  inquire  if  there  exists  a  utility  function  with  which  the  demand  function  is  associated.  If  yes, 
we  say  that  the  demand  function  is  integr&ble. 

We  will,  however,  restrict  the  utility  functions  considered  from  now  on  to  those  homogeneous 
in  v  of  degree  -f  1.  To  conform  with  this  definition,  we  redefine  the  utility  function  (1.1.2)  for  an 
individual  to  be  UJ(X)  =  —  [(SJ  —  X)'M3(S3  —  X)]1/3.  We  call  a  function  a  per  capita  utility 
function,  U (X)  =  —Z(v),  v  =  S  —  X,  if  (i)  Z(v)  is  twice  differentiable  for  all  v  (except  possibly 
at  v  =  0);  (ii)  Z[y)  is  a  homogeneous  function  in  v  of  degree  1  along  every  ray,  i.e.,  Z(av)  =  aZ(v) 
for  all  v,  a  >  0;  (iii)  Z(v)  is  strictly  convex  between  any  two  points  t;1  ^  v2  satisfying  the  budget 
constraint  p'v  =  1: 

A Z(v*)  +  pZ(v2)  >  Z( Av1  +  pv2)  for  all  (A  >  0,  p  >  0,  A  -I-  p  =  1)  ;  (2.2.0) 

and  (iv)  Z(v)  >  0  for  all  v  ^  0. 

We  first  show  that  Z(v)  under  the  second  definition  qualifies  as  a  utility  function  under  the  first 
definition.  We  then  seek  a  necessary  and  sufficient  condition  that  U(X)  exists  when  the  demand 
function  is  given  by  (2.1.0).  The  diagram  on  the  next  page  outlines  the  logical  dependence  of  the 
various  theorems  upon  one  another  which  yield  or  are  implied  by  this  condition. 
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Theorem  2.9 


Sufficiency,  Theorem  2.8 


Theorem  2.1.  min  Z(v)  subject  to  p'v  =  1  is  attained  at  a  Gnite  point  v  =  v*  and  is  unique. 
Denoting  v(p,  1)  =  v* ,  the  expected  per  capita  demand  function  associated  with  Z(v)  is  S  -  X  = 
v(p,  J)  =  J  •  v(p,  1)  where  J  =  (pS  -  Z). 

Proof.  Assume  on  the  contrary  that  inf  Z(v)  is  not  attained  at  any  finite  point.  In  this  case,  there 
exists  a  strictly  decreasing  sequence 

^(v1)  >  Z(v2)  >■•>  Z(v *)  >  •••  >  inf  Z{v)  (2.3.0) 

such  that  Z{y*)  — »  inf  Z(v)  >  0  as  J|v‘j|  — *  +oo,  and  such  that  p'v*  =  1  for  all  t.  A  subsequence 
can  be  chosen,  so  that  normalized  vectors  v*/||v*jj  — ►  v°,  where  v°  ^  0  and  Z(v°)  >  0.  For  this 
subsequence,  due  to  homogeneity  of  Z(y) , 

*(«*)  =  Ml  ■  2(»7II»‘II)  -  IM|  •  Z(.°)  -  +oo  ,  (2  3.1) 


contradicting  that  Z(v‘)  is  a  strictly  decreasing  subsequence. 


Therefore  Z(v)  attains  its  minimum  at  either  a  unique  finite  point  v*  or  attains  it  at  least 
two  finite  points  v1  /  t>2.  Now  v1  and  v2  cannot  be  on  the  same  ray  and  both  satisfy  the  budget 
constraint  p'v  =  1.  If  not  on  the  same  ray,  then  by  (2.2.0),  any  convex  combination  of  Av1  -f  pv2 
satisfies  the  budget  constraint  and  yields  a  lower  Z(v),  contradicting  Z(yx)  =  Z(v2)  are  minimum 
values. 

To  prove  v(p,  J)  =  J-v(p,  1),  following  the  same  proof  as  above,  we  know  that  min  Z(v)  subject 
to  a  budget  constraint  p'v  =  J  or  [p'/J)v  =  1  is  attained  at  a  unique  point,  v*  =  v(p,J).  Let 
v1  =  v(p,  1)  and  assume  on  the  contrary: 

Z(J  ■  v1)  =  JZ{yx)  >  Z(v*)  =  JZ{v*/J)  ,  (2.3.2) 

implying  that  v*/J,  which  satisfies  pv  =  1,  maps  into  Z(v*/J)  <  Z(vx)  a  contradiction.  | 

We  defined  pt  as  fixed  period  t  prices.  Let  p  be  a  rescaled  price  vector  that  maps  into  v  by  the 
relation:  v  =  G(p).  We  have  the  following  definitions  and  relations: 

S  -  X  =  u  =  G(p)  ,  (2.4.0) 

G(p)  =  Siip'IPp)-1  ■  H*p  (2.4.1) 

H(p)  =  £t(p‘  H'p)~x  ■  W  (2.4.2) 

p’v  =  1  for  all  p  ^  0  .  (2.4.3) 

Two  lemmas  about  G(p)  and  H(p)  : 

Lemma  2.1.  Elements  (k,  l)  of  the  matrix  H(p)  =  £i(p’Htp)~1  -H*  are  homogeneous  functions  of  p 
of  degree  -2;  H(p)  is  square,  symmetric,  positive  definite,  moreover  pf[H(p)\p  =  p'G[p)  —  p'v  =  1 
for  all  p  ^  0. 

Proof.  By  definition,  element  ( k,t )  of  H[p)  is  £i{jpf H% p)~x  H'kl  and  is  equal  to  element  {l,  k) 
since  H'kt  =  H^.  It  is  positive  definite  because  the  assumed  positive  definiteness  of  H'  im¬ 
plies  q‘[H(p)]q  =  £i(p' Hl p)-1  {q' H% q)  >  0  for  all  p  ^  0,  q  ^  0.  In  particular  p'\H(p)]p  = 
£<(p'H<p)-1(p,H‘p)  =  1.  I 

Lemma  2.2.  The  elements  of  the  vector  v  =  G(p)  =  H(p)p  =  £i(p'  H*p)~x  •  Hlp  are  homogeneous 
functions  of  degree  —1.  The  matrix  dv/dp'  is  square  and  symmetric. 

Proof.  We  define  element  (k,  t)  of  dv/djf  as  dvkj  dpt  where  vk  is  fc-th  component  of  t;  and  pt  the 
f-th  component  of  p.  From  (2.4.0),  (2.4.1): 

=  f.(p'^p)-1  •  Hip  ,  (2.5.0) 

where  Hxk.  denotes  the  Ar-th  row  of  Hl.  Clearly  vk  is  a  homogeneous  function  in  p  of  degree  —  1. 
Taking  partials, 

=  _  ( Hjp)(p'Hi  +  Hi,ph  dvi 

dpt  *  xp>Hip  ( p'fif’p )2  J  dpk 


(2.5.1) 


TO  TV\  RAW  1*1  jvi  *i\ 


where  the  f-th  column  H't  =  (Hlt.)'  and  H*kt  =  because  H'  is  symmetric.  | 


Theorem  2.2.  The  unique  v  that  minimizes  Z(v)  subject  to  the  budget  constraint  p'v  =  1  satisfies 
the  first-order  condition: 

i  ~  =  P  •  (2-6.0) 


Z  d v  r'  v  '  '  ' 

Conversely,  given  any  v  =  v* ,  there  exists  a  unique  p  such  that  min  Z  subject  to  p'v  =  1  is  attained 
at  u* ,  namely: 

4—  L...  =P-  (2.6.1) 


Z  ’  dv  P 


Proof.  .  Forming  the  Lagrangian  Z(v)  -  A  (p'v)  and  setting  its  partial  derivates  to  zero,  we  obtain 
dZ/dv  =  X  ■  p  where  A  is  chosen  so  that  p'v  =  1.  Applying  Euler's  Theorem  for  homogeneous 
functions  to  Z(v)  of  degree  +1: 

Z  =  v'(^)  =  v'(Xp)  =  A {v'p)  =  A  ,  (2.6.2) 

whence  (2.6.0).  Conversely,  given  v  =  v* ,  the  p  defined  by  (2.6.1)  satisfies  v'p  =  v'(dZ/dv)/Z  — 
Z/Z  =  1  by  Euler’s  Theorem.  Suppose  now  for  p  given  by  (2.6.1)  that  min  Z(v)  subject  to  p'v  =  1  is 
attained  at  some  other  v  =  v°  /  v*  satisfying  (2.6.0).  Then  since  Z(v)  is  convex  and  differentiable, 
both  n°  and  w*  satisfy  necessary  and  sufficient  conditions  (2.6.0)  to  be  global  minima  points, 
contradicting  uniqueness  established  in  Theorem  2.1.  | 


Theorem  2.3.  A  necessary  condition  for  the  existence  of  a  utility  function  U(X)  =  —Z(v),  asso¬ 
ciated  with  the  demand  function 

S  -  X  =  J  •  G(p)  ,  J  =  pS-I,  (2.7.0) 

is  that  the  inverse  function  of  v  =  G(p),  exists,  namely, 

p  =  G-1(v)  =  ~  ,  p'v-l  ,  v  =  S~X.  (2.7.1) 

Proof.  If  v  =  G(p),  then  by  (2.4.3)  p'v  =  1  for  all  p.  If  a  utility  function  U(X)  =  ~Z(v)  exists 
associated  with  the  demand  function  (2.7.0),  then  minimizing  Z(v)  subject  to  p'v  =  1  satisfies 


(2.6.0).  Since  Z  is  a  function  of  v,  (2.7.1)  states  that  p  can  be  expressed  as  a  function  of  t>  when 
p'v  =  1.  Hence  a  necessary  condition  for  the  existence  of  a  utility  function  is  that  v  =  G(p)  have 
an  inverse  function  p  =  G-1(v).  | 

Comment.  There  are  two  ways  that  the  function  v  =  G(p)  can  fail  to  have  an  inverse.  The  first 
is:  given  v,  there  exists  no  p  that  satisfies  the  equation  v  =  G{p).  The  second  way  is:  given  a 
particular  v,  there  is  more  than  one  p  satisfying  the  equation.  We  will  prove  that  first  way  can 


/V*.  ' 


w  v.  >r -  V'.'vjv.’w'f-Tfjvw.  rvsvjvjv.r^  •«■;  •. 


never  happen  and  the  second  way  can.  However,  under  certain  conditions  that  we  will  specify  later, 
there  is  a  unique  solution  for  all  choices  of  v  ^  0. 

Theorem  2.4.  Given  any  v  0,  there  are  always  one  or  more  p  ^  0  satisfying 

v  =  G(p)  =  &(p'iFp)_1  •  &P  ,  «#0.  (2.8.0) 


Proof.  Let  i  =  (1, . . . ,  n)  and  let  fl  =  (Aj  :  A*  >  0,  A*  =  1}  be  the  n  -  1  dimensional  simplex. 
We  now  define  a  continuous  mapping  A  — *•  p  €  fl:  Choose  any  A  e  ft;  let  if  =  A,#*;  determine 

p  =  nH~1v.  Note  that  H  is  positive  definite  because  it  is  a  convex  combination  of  positive  definite 
matrices  IP;  hence  H~x  exists;  hence  p  can  be  computed.  Next,  form  pi  =  (p'iPp)-1  >  0  for 
each  t  and  set  p<  =  pi/ Pk-  It  is  not  difiicult  to  show  that  the  mapping  A  — ♦  p  is  continuous. 
According  to  Brouwer’s  Theorem,  a  fixed  point  A0  — »  A0  exists.  For  this  A0  =  p  determine 


p  =  nH~1v  where  H  =  E"A°fP  , 

(2.8.1) 

Pi  =  A?  =  (p'Hip)~1/'Eli  (p'H’p)- 1  , 

(2.8.2) 

nv  =  Hp  =  E?A ?^p 

(2.8.3) 

=  [E?(p'H‘p)-1  •  /PpJ/EJfp'fTp)-1 

v  =  (l/njSJtp'^p)-1  •  /Pp  =  G(p)  . 

(2.8.4) 

where  p  =  p  •  E "(p'H’p)  1 ,  whence  (2.8.0).  | 

Theorem  2.5.  Given  certain  v  and  certain  IP,  it  is  possible  there  exist  more  than  one  p  satisfying 
v  =  G(p),  implying  G~1(v)  does  not  exist  in  general. 


Proof.  Let  n  =  2  where  $  =  1 . n  and  let  m  =  2  where  matrices  H*  are  m  x  m.  Let 


H1  = 


/  2.0 

\  -6 


-.e\ 
).2  /  ’ 


if3  = 


(  0.2 

V--6 


(2.9.0) 


In  the  case  of  n  =  2  and  m  =  2,  it  is  not  difiicult  to  show  that  there  are  three  different  p 
satisfying  v  =  (1/2)  (p'ff*p)-1  "ff’p  generated  by  three  different  choices  of  A°  and  corresponding 

Aj  =  1  -  A°,  see  (2.8.1)  and  (2.8.2).  One  of  these  three  values  of  A°  can  be  real  and  the  other  two 
conjugate  pairs,  a  ±  0\/-T,  or  all  three  can  be  real.  For  the  special  case  (2.9.0),  there  are  three 
real  solutions: 


p  =  (1, 1)  ,  p  =  (.48382,  1.51618)  ,  p  =  (1.51618,  .48382) 

Since  these  solutions  satisfy  v  =  G(p),  all  three  satisfy  p'v  =  1.  | 
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Theorem  2.6.  For  some  choices  of  H',  v  =  G(p)  has  an  inverse  p  =  G  1(v). 

Proof.  Assume  H*  =  R  for  all  i.  It  is  easy  to  verify  that  a  p  which  satisfies  (2.10.0)  is  given  by 
(2.10.1): 

v  =  (p'Hp)-1Hp .  (2.10.0) 

p  =  {v'R-'v)'1  •  H-'v  ,  (2.10.1) 

The  RHS  of  (2.10.1)  is  G~1( t/)  because,  as  we  will  now  prove,  this  solution  is  unique:  Assume,  on 
the  contrary,  p^p  also  satisfies  (2.10.0): 

v  =  (p' Rp)-1  Hp  .  (2.10.2) 

Substituting  this  expression  for  v  into  (2.10.1),  we  obtain  p  =  p,  a  contradiction.  | 

When  v  =  G(p)  has  an  inverse  p  =  G~1(v),  we  will  make  use  of  the  following  lemmas. 

Lemma  2.3.  v  =  G(p)  and  p  =  G-1(u)  are  homogeneous  functions  of  degrees  p  =  -1  and 
1/p  =  -1  in  p  and  v  respectively. 

Lemma  2.4.  If  the  matrix  dv/dp'  is  non-singular  at  p,  then  its  inverse  dp/dv'  exists. 

Proof,  dv/dv  =  identity  =  \dv/dp'\  •  \dp/dv')  from  the  theory  of  implicit  functions.  | 

Lemma  2.5.  The  matrix  dp/dv'  is  symmetric. 

Proof.  This  follows  from  Lemmas  2.2  and  2.4.  | 

Lemma  2.6. 

v  =  ~[dv/dp']p  •  v'  =  -p'[dv/5p']  ;  (2.11.0) 

p  =  -(dp/9v']v  ;  p'  =  -v'\dp/dv'\  .  (2.11.1) 

Proof.  The  first  part  of  the  Lemma  follows  from  Euler’s  Theorem  for  homogeneous  forms  of  degree 
—  1,  see  Lemma  2.3,  and  the  second  part  follows  from  the  first  part  and  symmetry,  see  Lemma  2.5. 

I 

Lemma  2.7. 

P'§pp+(p'v)2  =  0  (2-12.0) 

Proof.  By  (2.11.0),  the  first  term  is  -j/v  and  j/v  =  1  by  (2.4.3).  | 
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g'— ?  +  ( q'v )2  >0  for  all  q  ±  ap,  p  ^  0  . 
op 


Proof.  By  definition  of  a  utility  function,  its  correpsonding  disutility  function  satisfies  Z(v)  >0 
for  all  v  ±  0  and  d2Z/dv2  is  a  positive  semi-definite  matrix  such  that: 


,d2Z  n  ,d2Z  n  .  „  , 

v  g~2v  =  0  ,  u  -^-u  >  o  ,  for  all  u  ±  aw  . 


We  have  shown  by  (2.7.1): 


p  =  {dZ/dv)/Z 


(2.12.2) 


(2.12.3) 


At  any  v  ^  0,  dp/dv'  exists  since  dp/dv '  =  {\/Z)d2  Z/dv2  —  (1  /Z)2(dZ/dv)2  and 
Z(v)  ^  0.  From  d2Z/dv2  =  d[dZ/dv]/d v'  and  (2.12.3): 


^l  =  ^PZ-hp^  =  ^.z  +  (pP').z 

dv2  dv<  ^Pdv>  dv1 

=  M|£+W>'] 


(2.12.4) 

(2.12.5) 


By  (2.5.1),  the  inverse  of  dp/dv',  namely  dv/dp',  also  exists,  at  any  point  »^0,  p  /  0.  Therefore, 
for  u  ^  a»,  we  can  find  a  g  satisfying 


u=aPq' 

Z(v),  by  assumption,  is  strictly  convex  except  along  rays.  We  therefore  require 
°  <  'i'~dv*U  =  Z[U'fri?U  +  for  811  u  ^  av>  v  #  0  > 

a 

=  ^[?'^7  +  (7'v)3]  for  all  g  ^  a  •  p,  p^O, 
by  (2.12.6).  Whence  (2.12.1)  since  Z  >  0.  | 


(2.12.6) 


(2.12.7) 

(2.12.8) 


Theorem  2.8.  A  sufficient  condition  that  the  inverse  function  p  ~  G  1(v)  exists  when 
v  =  G(p)  =  £,  (p/f*p)-1  •  Hlp  is 


q'~q  +  {q'v)2  >  0  for  all  g  ,  g^a-p,  p^0. 
dp 


Proof.  We  first  show  that  the  function  W[p),  defined  by 


(2.13.0) 


log  W(p)  =  i  £  log(p'ff*p)  ,  p  ^  0  , 


(2.13.1) 


»»*  h"  i' i' »'  i ‘  ■,r  »' r l S**’  »»»!i*iJ  i 
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under  the  hypothesis  (2.13.0),  is  (a)  convex;  (b)  W(p)  >  0  for  p  ^  0;  (c)  its  Hessian  d2W /dp2 , 
is  strictly  positive  definite  for  all  directions  q  ^  a  •  p;  and  finally,  (d)  p'(d2W /dp2)p  =  0.  Clearly 
IV (p)  >  0.  The  first  partial  is: 


Letting  v  =  G(p),  dW /dp  =  Wv.  Applying  (2.13.3),  the  second  partial  is: 

wmW*&+ 

a2rir  <9*i 


(2.13.2) 

(2.13.3) 


(2.13.4) 

(2.13.5) 


by  hypothesis  (2.13.0)  and  IV  >  0.  Finally,  if  we  replace  q  by  p  in  (2.13.5),  then  by  Euler’s  Theorem, 
p'(dv/dp')p  =  p'(— v)  =  —1  and  (p'v)2  =  +1,  whence 

d2W 

_lu  rr _ a  { one! 


P  dp2  P  0  ' 


(2.13.6) 


Conditions  (2.13.5)  and  (2.13.6),  and  IV  >  0  imply  that  d2W /dp2  given  by  (2.13.4)  is  positive 
semi-definite. 

We  now  use  the  properties  of  W  (p)  to  show  for  any  given  v  there  exists  a  unique  p  such  that 
v  =  G(p).  Suppose  not  true  and  for  some  given  v  ^  0,  there  exists  two  vectors  p  —  p  and  p  =  p 
such  that  v  —  G(p)  is  satisfied.  Consider  now  the  problem  of 


min  IV  (p)  subject  to  v'p  =  1  , 


(2.14.0) 


where  v  is  given  and  p  is  variable.  Since  W(p)  is  convex,  we  find  its  minimum  at  some  p  =  p*  by 
setting  the  partials  of  the  Lagrangian  W  (p)  -  A(v'p)  to  zero  obtaining 


dW 

=  \v  =  Wv, 
dp 


(2.14.1) 


where  A  =  IV  because  W  =  p'(dW/dp)  —  p'( Av)  =  A  (p'v)  =  A.  Equation  (2.13.3)  and  (2.13.1) 
together  expresses  dW/dp  in  terms  of  p  while  (2.14.1)  places  a  condition  on  p  to  qualify  as  a 
minimizing  p*.  Equating  (2.14.1)  with  (2.13.3),  we  see  that  p*  satisfies  v  =  £,(p' H'p)~1  •  H'p,  and 
therefore  p*  =  p  and  p*  =  p  are  both  optimal. 

But  minimizing  of  a  convex  function  W  (p)  subject  to  linear  constraint  v'p  —  1  can  only  have 
global  minima.  Therefore,  W (p)  and  W (p)  are  both  global  minima  and  IV (p)  =  W (p).  Vectors  p 
and  p  cannot  lie  on  same  ray  and  both  satisfy  the  budget  constraint  (2.14.0).  Therefore,  by  (2.13.5), 
in  the  direction  q  =  p  -  p  from  p  to  p,  (p  -  py\d2W /dp\(p  -  p)  >  0  for  all  p  =  Xp  +  pp,  A  +  p  = 
1,  (A ,p)  >  0.  Because  of  strict  convexity  along  this  segment,  W(p)  <  W(p)  =  W (p)  holds;  also 
p'v  =  p'v  =  p'v  =  1  holds,  contradicting  p  and  p  both  being  minimizing  points.  | 
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Theorem  2.9.  A  necessary  and  sufficient  condition  that  there  exists  a  utility  function 
U (X)  =  —Z(v),  v  =  S  —  X,  associated  with  the  expected  per  capita  demand  function  S  —  X  = 
(p’S  -  I)£i(p’ H'p)*1  •  H'p  for  period  t  is 


d\) 

q  dp;q  +  (q>V^  >  °  ’  f°T  011  q  ^  “  * P  ’  P  ^  °  ‘ 

The  function  Z(v)  is  the  mapping  v  —*  Z  defined  by  the  following  procedure: 


(2.15.0) 


Step  1.  Find  the  unique  p  satisfying 


v  =  Siip'IPp)-1  ■  IPp 


Step  2.  For  any  fixed  scalar  7  >  0,  find  Z  satisfying 


log lZ  =  -\£i  log (p'iTp),  7  >  0  , 


(2.15.1) 


(2.15.2) 


Proof.  Necessity  of  (2.15.0)  has  already  been  demonstrated.  Therefore  we  assume  (2.15.0)  is  true; 
by  Theorem  2.8,  p  =  G-1(v)  exists,  hence  p  of  (2.15.1)  is  unique.  Our  objective  is  to  show  that 
Z(v)  =  'jZ(v),  7  >  0,  is  a  disutility  function  associated  with  v  =  G(p).  To  qualify  by  our  definition 
as  a  disutility  function,  we  must  show  properties  A,  B,  C,  D,  E  below: 

A.  Z(v)  >  0  for  all  t>  /  0.  This  follows  by  rewriting  (2.15.2), 


i  =  \(p'Hxp)  ■  (p'H2p)  ■ .  •  ( p'Hnp)}-1'2n  >  0  ,  for  all  p  /  0  . 


(2.15.3) 


B.  Z(v)  is  homogeneous  in  v  of  degree  =  1.  Proof:  Note  v  is  homogeneous  in  p  of  degree  —1  by 
(2.15.1)  and  Z  is  homogeneous  in  p  of  degree  -1  by  (2.15.3),  implying  Z  is  homogeneous  in  v  of 
degree  1. 

C.  Minimizing  Z(v )  subject  to  the  budget  p'v  =  1  implies  the  demand  function  v  =  G(p).  Proof: 
According  to  Theorem  2.3,  the  inverse  demand  function  associated  with  minimizing  Z(v)  subject 
to  p'v  —  1  is 

1  dZ  . 

<2I54> 

where  the  left  hand  side  is  viewed  as  a  function  of  v.  We  must  therefore  verify  in  the  case  when 
v  =  G(p)  and  Z  is  defined  by  (2.15.1)  find  (2.15.2)  that  the  left  hand  side  of  (2.15.4)  is  indeed 
G_1(v).  From  (2.15.2): 


±dz  =  a  log  z  =  (-i/2)a[£iog(p'jPp)]  dp 

Z  dv'  dv'  dpt  dv' 

=  -iwffvr1  •  n-p]'  .  ||£] 


(2.15.5) 


(2.15.6) 


s 


where  dp/dv'  refers  to  the  functional  relation  between  p  and  v  implied  by  v  =  G(p)  and  where 
the  last  step  follows  from  Euler’s  Theorem  for  homogeneous  forms,  see  (2.11.0).  Therefore,  the  p 
generated  by  (2.15.4)  is  the  same  one  that  is  related  to  v  by  v  =  G(p). 

D.  The  function  Z(v)  is  convex.  Proof:  Following  the  steps  of  Theorem  2.7, 


,d2Z  A.  ,  dv  .  .  dv 

u'fo2u=Z{q  dPq+(qu)  1  >0,  U=3^?’  *^ap>  p^°- 


(2.15.7) 


According  to  our  hypothesis  (2.15.0)  and  Z  >  0  by  (2.15.3),  the  Hessian  [d2  Z  /  dv2\  is  strictly 
positive  definite  in  the  direction  of  any  two  points  not  on  the  same  ray.  Note  (2.12.0)  holds  along 
a  ray. 

E.  General  Z(y)  =  lZ(v)  where  7  >  0  is  a  scalar  constant.  Proof:  By  Theorem  2.3  and  (2.15.4), 
both  general  Z[y)  and  Z(v)  satisfy 


Integrating 


J_dZ  _  1  dZ  _ 

Z  dv  ~  Z  dv  ~P 


log  Z  —  log  Z  "I-  log  7  ,  7  >  0 


(2.15.8) 


(2.15.9) 


whence  Z  =  7  •  Z.  From  now  on  we  will  assume  7  =  1  and  Z  =  Z.  | 


Theorem  2.10.  A  necessary  and  sufficient  condition  that  there  exists  a  utility  function  U[X)  = 
-Z(v),  v  =  S  -  X,  associated  with  the  expected  per  capita  demand  function  S  -  X  =  (p'S  -  T)  ■ 
<f,  (p'/f‘p)_1  •  H'p  for  period  t  is 


c.\£Ei  _  +  if.s^Ei*  >  0 

t'p/H'p  (p'H'pp'  ‘‘"t/H'p 


for  all  p  0,  q  ^  ap. 


Proof.  Substituting  into  (2.15.0),  v  =  f,(p'if*p)  1  •  H'p.  Noting  (2.5.1), 


dv  _  H* 

dp>  (p'H'p)2  1 


form  which  we  obtain  (2.16.0).  | 


(2.16.0) 


(2.16.1) 


PART  m: 

TWO  SUFFICIENT  CONDITIONS  FOR  INTEGRABILITY 


Having  found  a  necessary  and  sufficient  condition  that  a  per  capita  utility  function  exists, 
namely 


rcfWq  W&p)2.  r 


>  0 


(3.1.0) 


holds  for  all  p  ^  0,  q  a  •  p,  we  now  seek  conditions  on  H'  in  the  population  urn  that  guarantee 
this.  We  prove  two  important  theorems  which  show  that  the  distribution  of  H'  in  the  population 
urn  would  have  to  be  highly  polarized  as  in  the  example  of  Theorem  2.5  in  order  for  condition 
(3.1.0)  to  fail.  It  is  not  difficult  to  show  if  Hi  =  H  for  all  *  that  (3.1.0)  holds.  Let  H  be  any 
positive  definite  matrix,  for  example  H  =  CiH1  where  =  1.  One  measure  of  how  much  H' 

differs  from  H  is  to  form  H'  =  H~1!2 H' H~1!2  and  compare  H*  to  J,  the  identity  matrix.  The 
eigenvalues  of  I  are  all  unity  and  the  ratio  of  its  hightest  to  lowest  eigenvalues  is  unity.  Therefore, 
we  can  study  how  much  the  ratio  of  the  highest  to  lowest  eigenvalues  of  Hl  has  to  differ  from  unity 
before  condition  (3.1.0)  is  violated. 


Indeed  we  will  show  in  Theorem  3.4  that  this  ratio  would  have  to  be  greater  than  3  +  \/8  for 
each  »  in  order  for  condition  (3.1.0)  to  fail.  But  even  this  will  not  cause  failure  when  the  axes 
of  various  ellipsoids  pH'p  =  constant  are  randomly  rotated  to  some  extent  with  respect  to  each 
other;  if  the  rotations  are  uniformly  distributed  (or  nearly  so)  in  Rn,  then  we  will  show  in  Theorem 
3.5  that  regardless  of  what  the  ratio  of  highest  to  lowest  eigenvalues  of  H ’  are,  condition  (3.1.0) 
will  hold.  In  other  words  the  distribution  of  the  price-  cross  effect  matrices  H'  would  have  to  be 
exceptionally  highly  skewed  for  a  per  capita  utility  function  to  fail  to  exist. 


Theorem  3.1.  Let  H  be  any  positive  definite  matrix  and  let  Hl  ~  H~1!2  H'  H~1f2 ,  then  the 
necessary  and  sufficient  condition  that  a  utility  function  exists  is  equivalent  to  finding  conditions 
on  H%  so  that  for  all  p  ^  0  and  q  ap: 


^_2m]  +  [,t^£l2>  o 


V#’P  (p'H'p)2 


P*  H'p 


(3.1.1) 


Proof.  The  matrix  H1!2  is  not  unique.  There  is  a  way  to  choose  it  so  that  H 1/f2  is  symmetric, 
namely  Hx!2  =  ED1/2 E'  where  E  is  the  matrix  of  eigenvectors  of  H  and  D  is  the  diagonal  matrix 
whose  diagonal  is  the  eigenvalues  of  H.  For  properties  of  E  see  proof  of  Theorem  3.3.  Substituting 
p  =  H~1t2p,  q  =  H~l/2q  into  (3.1.0)  and  then  relabeling  ( p,q )  as  ( p,q )  in  order  not  to  have 
proliferation  of  symbols  we  obtain  (3.1.1).  | 

Comment.  If  H'  are  all  close  to  H  —  CiH*,  then  Hx  will  be  close  to  £</?*  =  7,  the  identity.  If 
H*  =  I  for  all  i,  condition  (3.1.1)  reduces  to  showing  p2q2  -  ( p’q )2  >  0  for  all  p  ^  0,  q  ^  0 ,p  ^  aq 
where  p2,q2  denotes  p'p,  q'p  respectively.  But  the  latter  is  always  true  because  it  is  the  same  as 
p2q2  sin2  9  where  9  is  the  angle  between  vectors  p  and  q. 


Theorem  3.2.  The  necessary  and  sufficient  condition  (3.1.1)  that  a  utility  function  exists  is  equiv¬ 
alent  to  finding  conditions  on  H '  so  that  for  all  p  ^  0,  g  ^  0,  p*g  =  0: 


rt=o. 


(3.1.2) 


Proof.  The  only  difference  between  (3.1.1)  and  (3.1.2)  is  the  requirement  that  g  be  orthogonal  to 
p.  Because  of  homogeneity  of  (3.1.1)  we  can  rescale  p  and  q  so  that  p2  =  1  and  g2  =  1.  Let  9  be  the 
angle  between  p  and  q  so  that  p'q  =  cos  9.  Note  that  the  sign  of  the  L.H.S.  of  (3.1.1)  is  the  same 
if  we  replace  q  by  —q  so  that  we  need  only  consider  0  <  0  <  x.  Therefore,  the  condition  q  /  ap 
translates  into  sin#  ^  0.  We  replace  variables  (p, q)  in  (3.1.1)  by  p, q,  and  9  where 

p  =  p,  g  =  gain# +  pcos#,  cos 9=  p'q  where  sin 0,40.  (3.1.3) 

It  is  easy  to  prove  when  p2  =  1 ,  g2  =  1  that 

ft  ?  =  0,  g2  =  1,  p2  =  l.  (3.1.4) 

Substitute  (3.1.3)  into  (3.1.1).  After  much  cancellation  of  terms  and  factoring  out  of  the  common 
factor  sin2  9  >  0,  we  obtain  (3.1.2).  Because  (3.1.2)  is  homogeneous  in  p  find  q,  we  no  longer  require 
p2  =  1,  g2  =  1.  | 


Comment.  The  last  bracket  expression  of  (3.1.2)  is  not  likely  to  contribute  much  to  the  positivity 
of  the  L.H.S  For  example,  if  H  =  CiH1  and  all  Hl  =  £iHi  =  I,  the  identity,  the  last  term  would 
vanish  because  p*g  =  0.  Therefore,  if  we  drop  the  second  bracket,  it  is  sufficient  to  only  consider 
conditions  on  Hl  that  guarantee  that  the  first  bracket  expression  is  positive.  As  an  extremum  or 
worst  case  scenario,  we  will  look  for  conditions  on  H *  that  will  guarantee  for  every  i  that  ,  the 
corresponding  pair  of  ith  terms  in  the  first  bracket,  is  positive: 


paw  >0 


for  all  p2  ^  0,  g2  7^  0  L,ich  that  p/g  =  0  .  (3.1.5) 


To  simplify  the  discussion,  we  assume  edl  the  eigenvalues  of  ffl  are  distinct.  If  not  they  could  be 
made  so  by  a  slight  perturbation  of  Hx. 


Theorem  3.3,  a  sufficient  condition.  A  utility  function  exists  if  for  each  i: 


_  q'Dq 

p'Dp 


Aq’Dp )2 

(p'Dp)2 


>  0 


for  all  p'q  =  0,  p2  =  1,  g2  =  1 


(3.2.0) 


where  D  =  D'  is  the  diagonal  matrix  whose  diagonal  elements  are  <  ■  ■  ■  <  ,  the  m  distinct 

eigenvalues  of  B*. 


Proof.  All  the  eigenvalues  of  the  positive-definite  matrix  S *  are  positive.  Let  these  be  d*k  >  0. 
To  simplify  the  notation  we  will  write  dk  —  d^  and  denote  their  corresponding  eigenvectors  by 
Ek  -  E*k  rescaled  so  that  {Ek)2  =  1.  Then  by  definition 

HiEk  =  dk-I -Ek\  E‘lHi  =  dt’E'rJ  ,  (k,l)  =  1, ...  ,m  (3.2.1) 

where  7  denotes  the  identity  matrix.  Therefore  Ftff’Ek  =  dk  •  Ek  =  df  E?t  •  Ek-  Since  dk  ±  dt 
for  all  ly£k,  it  follows  that  E,tEk  =  0.  Therefore  the  matrix  of  eigenvectors  E  =  (f?i,  E^, . . . ,  Em) 
is  an  orthoncunal  matrix,  i.e.,  E'E  =  I,  and 

E'HiE=D.  (3.2.2) 


Substituting  q  =  Eq,  p  =  Ep  into  (3. I.E)  and  rescaling  so  that  p2  =  1,  q2  —  1,  p'q  =  0  we 
obtain  (3.2.0).  | 

Theorem  3.-*.  A  utility  function  exists  if  *,  r  every  i  the  diagonal  matrix  D  =  D'  of  eigenvalues 
of  H'  hw°  the  property  mint,-  >  0  where 

'•=  1-2^T55>°.  t3-3.0) 

where  pi  is  the  ratio  of  the  highest  to  lowest  eigenvalues  of  H',  or  equivalently  when 

Pi  <  3  +  2y/2  =  5.83  (  3.1) 


Proof.  Condition  (3.3.1)  follows  by  rewriting  (3.3.0)  as  (-p2  +  6 p  -  l)/(p  +  l)2  and  determining 
the  range  of  p  where  y  =  —p2  +6p  -  1  >  0.  To  prove  (3.3.0),  we  minimize  ti  given  by  (3.2.0)  subject 
to  p3  =  1,  q2  =  1,  pfq  =  0.  For  the  p,  q  that  yields  this  minimum,  let 

p' Dp  =  a  >  0  and  q'Dq  —  f)  >  0  ,  p'q  =  0  .  (3.3.2) 

Under  the  conditions  p' Dp  =  o  and  q'Dq  =  minimizing  U  reduces  to 

ma x(q'Dp)2  subject  to  (3.3.2)  .  (3.3.3) 

From  the  properties  at  the  maximum,  we  will  derive  a  relationship  between  ft  Dp  and  p2,  and  q'  Dq 
and  q2  that  will  allow  us  to  minimize  U  subject  to  j/q  =  0,  p2  =  1,  q2  =  1,  see  (3.2.0). 

Fact,  q' Dp  ^  0  at  the  maximum.  Proof:  We  assumed  diagonal  D  satisfies  0  <  d2  <  ■  ■  •  <  dm. 
Given  a  >  0  and  /?  >  0,  it  is  easy  to  find  p  and  q  satisfying  (3.3.2)  and  q' Dp  ^  0.  Therefore,  because 
everything  is  bounded  and  continuous,  a  maximum  exists  with  (q' Dp)2  >  0.  We  can  now  form  the 
Lagrangian 


L  =  (q1  Dp)2  -  X(p' Dp)  -  p{q'Dq)  -  2u(p'q) 


(3.3.4) 


and  set  dL/dp  —  0  and  dL/dq  =  0,  obtaining 


{q'Dp)  •  Dq  =  \  •  Dp  +  t/q  , 
{q'Dp)  -  Dp  =  p‘  Dp  +  vp  , 


(3.3.5) 

(3.3.6) 


Fact:  v  jL  0.  Proof:  Since  q'Dp  £  0,  it  follows  v  ^  0  because  v  —  0  would  imply  by  (3.3.5) 
that  p  and  q  are  proportional  and  p'q  ^  0,  contrary  to  hypothesis. 

Fact:  A  >  0,  p  >  0.  This  can  be  seen  by  multiplying  (3.3.5)  and  (3.3.6)  by  p  and  q  respectively  : 


{q'Dp)2  =  \a  =  pP  . 

Multiplying  (3.3.5)  and  (3.3.6)  by  q  and  p  respectively: 

{q'Dp)p  =  \{q' Dp)  +  t/q2  =  (q'Dp)3 /a  +  uq2 
{q'Dp) a  =  p{q'Dp)  +  vp2  =  {q'Dpf /P  +  vp2  . 


Hence 


{q'Dp)ctP  =  {q'Dpf  +  <xvq2 
{q'Dp)ap  =  {q'Dpf  +  Pup2 


implying 


(pa/a)  =  {q2/p)  and  a  —  P  when  p2  =  g2  =  1 


.2  _  „2 


We  now  rewrite  (3.3.5)  and  (3.3.6) 


[{tfDp)  •  D  -  ul]q  =  A  Dp 
\{q'Dp)  ■  D  -  vl)p  =  pDq 


(3.3.7) 


(3.3.8) 

(3.3.9) 


(3.4.0) 

(3.4.1) 

(3.4.2) 


(3.4.3) 

(3.4.4) 


where  I  is  the  identity  matrix.  Solving  (3.3.7)  for  A  and  (3.3.8)  for  u  and  substituting  into  (3.4.3): 


{(,'Dp)  D--~  [ap{q'Dp)  -  {q'Dpf}  • I}q=  •  Dp  , 


(3.4.5) 


and  an  analogous  expression  if  we  solve  (3.3.7)  for  p  and  (3.3.9)  for  v  and  substitute  into  (3.4.4). 
We  can  now  factor  out  {q'Dp)  ±  0,  obtaining  (3.4.6),  and  by  analogy  (3.4.7): 


{aq2D  -  \ap  -  {q'Dpf]  I}q  =  q2{q'Dp)  •  Dp 
{Pp* D  -  [ap  -  {q'Dpf  ]l}p  =  p2{q'Dp)  •  Dq 


(3.4.6) 

(3.4.7) 


where  the  expression  in  brackets,  are  the  same  in  (3.4.6)  and  (3.4.7)  since  aq2  =  Pp2.  Noting  that 
the  product  of  diagonal  matrices  is  commutative,  multiply  (3.4.7)  by  q2{q'Dp)D  and  interchange 


$ 

ft 

K*. 

& 


$ 

I 


tivii 


I 

& 

| 

■ 


on  the  L.H.S.  the  order  of  q2(q'Dp)D  and  the  bracket  term;  finally  substituting  L.H.S.  of  (3.4.6) 
for  q2(q' Dp )  •  Dp,  we  obtain 


{aq7D  -  [aft  -  (q' Dp)2]I}2q  =  p2q2(q'Dp)2  D2q  . 


Expanding  and  rearranging  terms  and  noting  ( aq 2)2  =  plq2aft: 


(3.5.0) 


\p2q2{a()  -  (q'Dpf)D2  -  2(aq2){af)  -  (jDp)2}D  -  {a/3  -  ( q' Dp)2)2  J\q  =  0  .  (3.5.1) 

We  observe  that  the  common  factor  {a/9  -  (q‘  Dp)2}  ^  0  by  (3.4.0)  or  (3.4.1)  and  the  fact  that 
v  ^  0.  Therefore  factoring  out  (aft  -  (q' Dp)2}  ^  0  from  (3.5.1),  we  obtain  (3.5.2)  and  an  analogous 
expression,  (3.5.3)  after  noting  aq 2  =  ftp2: 

\p2q2D2  -  2 (aq2)D  -  (aft  -  ( qDp)2)I\q  =  0  ,  (3.5.2) 

[p2q2D2  -  2 (aq2)D  -  (aft  -  (qDp)2)I\p  =  0  .  (3.5.3) 


Vector  relation  (3.5.2)  holds  for  each  component  qr*  of  q.  Since  q2  ^  0,  one  or  more  components 
qk  0.  Let  qk  t^O  for  some  k,  then  for  this  k: 

(P3?2)^  “  2 (<*q2)dk  -  [aft  ~  (q'Dp)2}  =  0  .  (3.5.4) 

Let  y  =  (p2q2)x2  —  2 (aq2)x  —  [aft  -  (q1  Dp)2]  be  a  parabola  expressing  y  as  a  function  of  x.  Now 
y  =  0  can  only  hold  for  at  most  two  values  of  *  =  dk,  i.e.,  for  say  k  =  1  and  k  —  m.  And  from 
(3  jnly  the  same  k  =  1  and  k  =  m  can  possibly  have  pk  ^  0,  the  case  of  only  one  k  being 
rule  jut  because  say  k  =  1  only,  then  both  pi  /  0,  qi  0;  and  all  other  p,  =  0,  y,  =  0.  But  this 
contradicts  p'q  =  0.  Representing  these  two  components  of  p  and  q  in  polar  coordinates  and  noting 


p'q  =  0: 

Pi  =  +IIpIIcos®i  Pm  =  +||p||  sin0  ,  and  p,  =  0  for  1  <  »  <  m  ,  (3.6.0) 

9i  =  — 1|?|| B\n0,  qm  =  +||g||cos0  ,  and  =  0  for  1  <  «  <  m  ,  (3.6.1) 

subject  to  p'Dp  =  a,  q'Dq  =  ft,  p'q  =  0.  Therefore 

a  =  p'Dp  =  p3[di  cos2  9  +  dm  sin2  0]  (3.6.2) 

ft  =  q'Dq  =  q2[d\  sin2  0  +  dm  cos2  0]  .  (3.6.3) 

Whence,  noting  p2 /a  =  q2 /ft,  by  (3.4.2), 

di  cos2  0  +  dm  sin2  0  =  di  sin2  0  +  dm  cos2  0  ,  (3.6.4) 

(dm  -  di)[cos2  0  -  sin2  0]  =  0  ,  (3.6.5) 


2' 


and  therefore  since  dm  ^  d\ , 


s3*  =  *, 


!2fl=  1 


a  =  p7[di  +  dm]/2,  P  =  q2[di+dm]/ 2, 
{qDp)2  =  pV  sin2  0  cos3  fl(dm  -  di)2  = 

(«i  +  dm; 


Recalling  a  =  fi  from  (3.4.2), 


0  ap(dm-d^  f  Jdm-d1)2 
a  a2(dm  +  di)3  (dm  +  di)2 


(3.6.6) 

(3.6.7) 

(3.6.8) 


(3.6.9) 


We  left  open  the  question,  which  pair  of  the  distinct  eigenvalues  di  <  d2  <  •  •  •  <  dm  to  select;  we 
show  now  that  the  best  choice  for  mint,  is  in  fact  the  ones  we  chose,  namely  d\  and  dm.  Let  j  >  t 
and  y  =  (dj  ~  di)/(dj  +  d,),  then  dy/ddj  =  2d,/(dy  +  d<)2  >  0  and  dy/ddi  =  — 2dy/(dy  +  d,)2  <  0. 
Therefore  the  ratio  increases  as  dy  — »  dm  and  also  as  d,-  — >  d\ .  Hence  the  lowest  2,  is  obtained  for 
«  -  1,  and  j  =  m. 


From  (3.6.0)  and  (3.6.1),  optimal  px  =  cos(ir/4),  pm  =  sin(jr/4),  and  pi  =  0  for  1  <  t  <  m;  and 
also  q i  =  -  sin(*/4),  qm  =  cos(jt/4),  and  $  =  0  for  1  <  »  <  m.  Therefore  the  minimal  solution  to 
(3.2.0)  is 

*  =  [‘  -  2wr+!$]  (3'70) 

Let  p  =  dm/ di  be  the  ratio  of  the  highest  to  lowest  eigenvalue,  then 


1- 


(3.7.0) 


*.  =  1  -  2|^n)5  =  (V  +  6P  -  !)/(/» +  !)2 


which  is  positive  for  all  p  in  the  range 


(3  +  2v/2)-1  <p<3  +  2V/2  =  5.83.  | 


(3.7.1) 


(3.7.2) 


Comment.  In  order  for  a  utility  function  to  exist  it  is  sufficient  that  £;t,  >  0.  Therefore  the 
condition  that  mint,-  >  0  for  each  i  is  a  worse  case  scenario  that  is  far  too  stringent.  We  know  for 
case  of  a  population  of  two  people  that  it  is  possible  for  &U  <  0,  and  this  is  true  even  if  the  term 
[£i(q' H'p)/(p' H'p)]2  is  added  to  U,  see  counter  example  of  Theorem  2.5.  In  a  large  population 
such  as  U.S.  we  would  expect  the  set  of  representative  ellipsoids  p'lPp  =  constant  in  the  urn  which 
can  be  rotated  into  one  another  to  have  to  some  extent  randomly  distributed  orientations  of  their 
axes.  If  we  rotate  all  such  ellipsoids  to  Hx  for  the  lowest «,  say  H,  the  effect  is  to  rotate  the  p  and 
q  to  a  random  position.  Therefore,  we  need  to  consider  for  each  such  H ,  the  average  or  expected 
value  of 

/(*,)  =  £-"[^“2(£f;)!]  '  e’  =  l,  ?  =  p'i  =  o,  (3.7.3) 


WOT 


$ 

■4 
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WOTi 


where  p  and  q  are  “randomly”  distributed  in  some  way.  We,  of  course,  do  not  know  what  is  the 
true  distribution  of  orientations  of  the  axes  of  the  ellipsoids  j/H'p  =  constant  that  can  be  rotated 
into  one  another.  We  will  prove  a  theorem  when  these  orientations  are  uniformly  distributed  about 
the  origin  in  fJm.  While  this  assumption  of  “uniformly  distributed”  is  not  realistic,  the  purpose  of 
Theorem  3.5  is  to  illustrate  that  the  more  the  orientations  vary  relative  to  one  another,  the  higher 
the  bound  p  <  3  +  2>/2  +  A  can  be  where  A  — »  +oo. 

Theorem  3.5.  If  the  orientations  of  the  axes  of  the  ellipsoids  p,Htp  —  constant  that  can  be 
rotated  into  any  particular  H*  =  H  are  uniformly  distributed  (or  nearly  so),  then  £Plqf{p,q)  >  0 
where  / (p,  q)  is  defined  by  (3.7.3);  moreover  Sit,  >  0,  implying  that  the  expected  per  capita  utility 
function  exists. 


Proof.  The  proof  consists  in  partitioning  the  set 


S(p,q)  =  {p,?|p3  =  1,  «3  =  1,  p'q  =  0} 


(3.8.0) 


of  all  admissible  (p,q)  into  equally  probable  subsets  T(R),  each  subset  obtained  from  T(I)  by  a 
general  rotation  R  in  Rm.  We  will  show  that  each  T(R)  has  the  property,  that  the  average  value 
of  f{p,q)  found  by  integration  over  T(R)  is  unity.  It  then  follows  that  the  integration  of  f(p,q) 
over  S(p,q)  will  also  be  unity.  As  noted  one  of  the  subsets  T(R)  is  the  set  T[I)  which  is  defined 
to  be  the  set  of  all  p  =  (cos  <f>, sin  <f>, 0, . . .  ,0)'  and  q  =  (-  sin  <P, cos  <f>, 0, . . .  ,0)'  for  0  <  <f>  <  2x. 

This  “reduces”  the  subproblem  to  considering,  in  place  of  H,p,q,  the  truncated  matrix  and 
vectors: 

H  =  "?12  ^  p  =  (cos  <f>,  sin  <p) ,  q  =  (—sm<f>,cos<f>)  .  (3.8.1) 

12  **22  J 

Let  the  matrix  of  eigenvectors  of  H  be  E  so  that  E1  HE  =  6.  The  diagonal  matrix  whose  elements 
are  the  eigenvalues  of  H,  namely  6(1,1)  =  dj  >  0  and  6(2,2)  =  da  >  0.  The  matrix  E  is 

orthonormal,  hence  we  can  rotate  p  and  q  in  R2  space  by  p  =  Ep  and  q  =  Eq  and  «/>  =  0  +  a  so  that 

the  axes  of  the  ellipsoid  p'Hp  =  constant  are  parallel  to  the  coordinate  axes.  After  this  rotation, 
our  problem  reduces  to  showing  for  T(7): 

2  fK/2  r di  sin2  6  +  da  cos2  9  „(d2  -  dx)2  sina  6  cos3  fll  Jg  _  ^ 

*  Jo  Ldi  cos3  0  +  da  sin3  9  [di  cos3  0  +  da  sin3  9)2  . 

where  p  =  (cos 0,sin  tf), q  =  (- sin cos 0).  Note  we  claim  that  relation  (3.8.2)  holds  independent 
of  what  the  diagonal  elements  of  6  happen  to  be. 

We  will  prove  (3.8.2)  in  a  moment.  We  can  obtain  the  subset  T(R),  by  rotating  the  subset 
T(7)  by  a  general  rotation  f2-1  about  the  origin  in  Rm  space.  Conversely  we  can  rotate  T(k)  into 
T(7)  by  p  =  Rp,  q  =  Rq,  where  R  is  an  orthonormal  matrix.  The  effect  is  to  replace  H  by  R HR 
and  carry  out  the  integration  over  T(7)  with  respect  to  RHR,  instead  of  H,  again  obtaining  unity. 
Since  each  rotation  R  is  equally  probable,  this  means  £P,qf(p,  q)  =  1  for  (p,  q)  G  S(p,  q). 


I 


We  now  show  (3.8.2)  is  true.  Substituting  *  =  cos2  6  +  p  sin20,  and  letting  da  =  pdi,  and 
Q  =  -p  +  (p  +  1)*  -  x 2,  (3.8.2)  reduces  to  proving 


l  ry+i-*  o1/a 


2rp  +  l  fp  dx  1  r  dx  [pQll2dx 


-ft  ¥£ 


xQ 'I* 


1  f'  dx  fp  Ql<2dx-[ 

2  A  Q1/2  Ji  x2  J-1 


(3.8.3) 


Denoting  a  +  6x  +  c*2  —  Q  =  ~p  +  (p  +  l)x  -  z2  and  noting  (62  -  4 ac)1/2  =  p  -  1,  we  obtain  from 
a  table  of  standard  integrals: 


rp  dx  1  .  _x  /  6z  +  2o  \ 

A  ^QTTa  -  (-0)1/2  S1D  ^x(62  _  4ac)i/2 ) 


1  niIt-1^+1)a?~2^r  __1_ 

-^/2am  v  x(p-i)  ;L=1-pV3 

fp  dx  1  .  /  — 2cz  -  6  \ 

Ji  Qi72  =  (-c)i/2  sm  y  (&2  —  4qc)i/2  y 

=-(^)L=' 

r  Q1/2,  Q1/2 1'  b  [p  dx  fp  dx 

Ji  x 2  “  x  L1  +  2;1  xq i/2+cy1  qi/2 


(3.8.4) 


(3.8.5) 


(3.8.6) 


Substituting  these  evaluations  into  (3.8.3),  yields  £p,qf(p,q)  =  1,  implying  =  1  >  0,  which 
is  sufficient  for  the  expected  per  capita  utility  function  to  exist.  | 

In  this  part  we  explored  how  close  the  price  cross-effect  matrices  Hl  must  be  to  a  typical  H  to 
guarantee  integrability.  One  measure  of  closeness  is  the  ratio  pi  of  the  highest  to  lowest  eigenvalues 
of  H'  =  It  is  sufficient  if  p,  <  3  +  2%/2  for  all ».  Another  measure  is  how  spread 

out  is  the  distribution  of  orientations  of  the  axes  of  the  ellipsoids  p’H'p  =  constanU,  the  more 
evenly  spread  the  better.  We  found  that  the  bound  on  p<  is  higher  the  more  the  tilts  of  the  axes  of 
the  ellipsoids  p’H'p  =  conatanti  are  uniformly  distributed  relative  to  one  another.  On  the  other 
hand,  if  individuals  t  in  the  population  tend  to  be  highly  polarized  as  to  their  consumption  tastes, 
as  in  the  example  of  Theorem  2.5  with  p\  =  pa  =  91  a  per  capita  utility  function  may  not  exist. 


PART  IV: 

INTEGRALITY  OF  THE  MULTIPLE-PERIOD  EQUILIBRIUM  MODEL 

In  Part  III  we  established  that  the  expected  per  capita  for  period  t  over  the  income  range  of 
interest  is  integrable  when  the  ratio  p,  of  highest  to  lowest  eigenvalue  of  Hl  is  less  than  3  +  V8 
for  all  i,  or  less  than  a  higher  bound  for  p,  if  the  axes  of  ellipsoids  p'H'p  =  const,  are  randomly 
disposed  to  one  another.  From  now  on  we  assume  that  a  utility  function  for  each  period  exists,  and 
consequently  an  inverse  demand  function  exists  except  for  a  scale  factor  to  be  determined.  This 
does  not  mean,  however,  in  the  context  of  a  time-staged  model  for  the  whole  economy,  that  there 
exists  an  objective  function  which,  when  maximized  subject  to  physical  flow  constraints,  implies  the 
additional  equilibrium  conditions  involving  prices  that  must  hold  between  the  producer /investor 
and  final  consumer.  It  depends  on  how  the  numeraire  for  prices  in  the  rate-of-return  formulas  are 
defined.  It  turns  out  that  when  the  numeraire  for  normalizing  prices  is  suitably  defined,  a  concave 
objective  function  for  the  economy  exists. 

Aggregate  quantities  corresponding  to  per  capita  quantities  are  denoted  by  bold  face  letters 
Xt,St,L,Ut,  Zt.  They  are  obtained  from  Xt,St,  It,Ut>Zt  by  multiplying  by  Pt,  the  size  of  popu¬ 
lation  in  period  t. 

The  mathematical  formulation  of  the  multi-period  model  is  along  the  classical  lines  of  an 
Arrow-Debreu  [1,5]  or  Scarf  [17]  equilibrium  model,  with  no  surprises  except  perhaps  for  the  inter¬ 
pretation  of  the  profitability  constraints  of  investors  as  rate-of-return  formulae  for  selecting  among 
different  investment  possibilities.  A  typical  production/investment  activity  j  in  period  t  has  a 
column  consisting  of  three  sets  of  fixed  coefficients  [R<(j),  —At(j),  —Dt[j)]'  per  unit  level  of  activ¬ 
ity  where  Bt(j)  is  the  input/output  vector  of  capacities,  resources,  and  flows  of  all  items  (goods) 
needed  for  production  and  capacity  formation  in  period  t  except  final  consumer  items;  At(j)  is 
the  output/input  vector  of  fined  consumer  items  in  period  t ;  and  Dt (j)  is  the  output  vector  of 
capacities,  resources,  and  intermediate  goods  left  over  or  produced  in  period  t  for  period  t  +  1.  The 
model  is  then  defined  by  five  sets  of  relations  numbered  (4.1.0)  through  (4.5.0).  The  first  two  we 
refer  to  as  the  primal  or  physical  flow  constraints 

In  (4.1.0)  below,  Yt  >  0  is  the  vector  of  aggregate  production  and  investment  levels  to  be 
determined  in  period  t.  In  words,  (4.1.0)  states  that  the  capacity,  resources,  intermediate  items 
jBtYt,  required  for  production  and  capacity  formation  at  levels  Yt  cannot  exceed  the  amount  of 
these  items  left  over  or  produced  by  period  t  -  1  activities  for  period  t  plus  the  vector 

kt  of  these  items  exogenously  supplied: 


Bt Yt  <  Z?«_iYt_i  +  kt  ; 


corresp.  dual  <rt  >  0  . 


(4.1.0) 


The  vector  of  dual  prices  corresponding  to  (4.1.0)  is  denoted  by  ot  >  0,  and  the  slack  vector  which 
turns  (4.1.0)  into  an  equation  is  denoted  by  bt  >  0.  According  to  the  theory  of  Arrow-Debreu, 


.'.O 


production  levels  Yt  and  prices  at  will  adjust  until  at  equilibrium  the  complementary  slackness 
conditions  a'tct  =  0  hold. 

In  (4.2.0)  below,  Xt  >  0  is  the  vector  of  aggregate  final  consumption  in  period  t  measured  in 
physical  units.  In  words,  (4.2.0)  states  that  consumption  cannot  exceed  AtYt,  the  net  output  from 
production  after  investment  less,  ft,  any  fixed  demand  (like  government  plus  any  minimum  floor 
provided  to  the  final  consumer  and  not  paid  for  by  income  for  consumption): 

— A<Yt  -f-X<  <  -ft  ;  corresp.  dual  jr*  >  0  .  (4.2.0) 

The  corresponding  vector  of  normalized  discounted  dual  prices  is  denoted  by  irt  >  0,  and  the  slack 
vector  which  turns  (4.2.0)  into  an  equation  is  denoted  by  srt  >  0.  At  equilibrium  r(  =  0. 

Dual  Constraints:  We  assume  that  a  numeraire  has  been  selected  in  each  period  relative  to  which 
prices  for  various  goods  are  measured.  If  so  we  say  the  prices  it  are  “normalized”.  Prices  <rt  >  0 
and  >  0  for  the  dual  constraints  are  defined  to  be  discounted  period  t  prices;  moreover  irt  is 
defined  to  be  discounted  normalized  period  t  prices.  Thus  itt  =  5t_1  ft  where  the  discount  factor 
is  5t_1  and  rt  is  the  vector  of  normalized  prices.  The  vector  of  unnormalized  period  t  prices  of  final 
consumer  items  is  denoted  by  pt .  We  have  no  need  for  a  symbol  for  prices  on  capacity,  reserves, 
and  intermediate  goods  relative  to  unnormalized  prices  pt  but  prices  on  these  items  relative  to 
normalized  prices  itt  are  denoted  by  ot. 

Relation  (4.3.0)  below,  which  we  will  derive  in  the  comment,  states  that  investors  must  receive 
at  least  their  minimum  rate  of  return  r  =  8~l  —  1  or  they  won’t  invest: 

A{jrt  <  -D'tot+ j  ;  corresp.  primal  Yt  >  0  .  (4.3.0) 

The  slack  vector  which  turns  (4.3.0)  into  an  equation  is  denoted  by  Yt  >0.  It  is  in  1  to  1 
correspondence  with  primal  variables  Yt  >  0.  At  equilibrium  YJYt  =  0. 

Comment:  Typically,  the  rate-of-return  formula  for  a  producer/investor  j  in  period  t  is  formed 
by  “pricing-out”  the  input-output  vector  of  physical  flows  per  unit  level  of  investment  activity  j, 
namely  [Bt(j),  -Dt(j)]',  and  then  looking  for  a  discount  factor  8  so  that  the  discounted 

cash  flows  is  zero: 

-p-'&'tBtU)  +  S^AU)  +  Ptr'^Dtti)  =  0.  (4.3.1) 

Note  the  use  of  normalized  prices  xt  (and  ot  relative  to  r(),  and  not  unnormalized  prices  pt.  This 
makes  the  rate-of-return  formula  inflation  free  meaning  that  the  investor  decides  what  his  minimum 
rate  of  return  r  should  be  without  multiplying  it  by  some  factor  for  future  inflation.  It  is  assumed 
that  an  investment  must  have  a  rate  of  return  not  less  than  a  minimum  rate  r,  or  a  d  >unt  factor 
6  <  S  =  (1  +  r)-1.  This  implies  that  investors  will  not  consider  investment  opportunity  j  if  its 
8  >  S.  On  the  other  hand  if  its  6  <6,  the  investment  is  profitable  and  the  level  Yt(j)  will  tend  in 
the  “real  world”  to  increase  indefinitely.  This  puts  pressure  on  capacities  and  resources  that  are  in 


short  supply,  thus  forcing  prices  and  S  up,  until  at  “equilibrium”  S  1  =  1  +  r.  We  express  these 
conditions  by  defining  Y <(t)  >  0  by 

-S'-'fi'tBtti)  +  +  *‘<r;+1Dt(i)  +  *t(j)  =  0  (4.3.2) 

and  requiring  Y t(j)  •  Y t(j)  =  0.  Substituting  St~1at  =  at,  =  xt,  P&t+i  =  <ft+i  and 

dropping  the  slack  vector  Y(j)  >  0,  this  relation  is  the  same  as  (4.3.0).  End  of  comment. 

In  (4.4.0)  below,  Ft (Xt )  is  the  inverse  of  the  expected  aggregate  demand  function  which  we 
assume  now  exists  except  for  a  scale  factor  to  be  determined  that  generates  normalized  period  t 
prices  i*  required  by  the  rate-of-  return  formula  (4.3.0): 

—xt  +  6t_1  ft(X<)  <  0  ;  corresp.  primal  X<  >  0  .  (4.4.0) 

The  corresponding  set  of  primal  variables  is  Xt  >  0.  When  X<  >  0,  which  is  usually  the  case,  this 
relation  becomes  an  equation.  The  slack  vector  that  turns  (4.4.0)  into  an  equation  is  denoted  by 
Xt  >  0.  At  equilibrium  XJX<  =  0.  The  case  when  the  inverse  demand  function  does  not  exist  will 
be  commented  on  in  a  moment. 

The  Complementary  Slackness  Conditions  are  that  all  variables  be  nonnegative  and 

a[ot  =  o,  *{*t  =  o,  YjYt  =  o,  x;xt  =  o.  (4.5.0) 

The  model  has  t  =  1,2, . . .  ,T  periods.  For  period  1  =  1,  the  term  D0Y0  is  omitted  in  relation 
(4.1.0).  For  period  t  =  T,  the  term  omitted  in  relation  (4.3.0). 

This  completes  the  mathematical  statement  of  the  time-staged  model.  The  remainder  of  Part 
IV  is  concerned  with  deriving  the  functional  form  of  the  inverse  aggregate  demand  function  J Ft(X<), 
and  the  utility  function  U(Xi, . . .  ,X()  for  the  full  economy  when  it  exists.  Before  doing  so  we  note 
that  if  the  economy  were  driven  by  a  utility  function  of  the  form  U  =  E^'-^UifXt),  the  Kuhn- 
Tucker  conditions  derived  by  maximizing  U  subject  to  primal  physical-flow  conditions  (4.1.0)  and 
(4.2.0)  would  give  rise  to  conditions  (4.3.0),  (4.4.0)  and  (4.5.0)  where  Ft(Xt)  =  dUt/dXt,  see 
reference  [12].  If  the  latter  conditions  hold  for  all  t  we  say  the  vector  functions  F((Xt)  in  the 
context  of  the  full  model  are  integrable.  If  not,  we  say  the  model  is  non-  integrable  and  no  utility 
function  for  the  economy  exists. 

The  dynamic  equilibrium  problem  is  well  defined  even  for  the  case  that  the  inverse  demand 
function  for  each  period  does  not  exist.  In  place  of  (4.4.0),  we  could  state  (a)  the  direct  aggregate 
demand  function  (4.6.0)  below  replacing  I(  by  jt[X<,  and  (b)  the  condition  prices  must  satisfy  when 
they  are  normalized  and  discounted.  It  is  outside  the  scope  of  this  paper  to  discuss  whether  or  not 
an  equilibrium  solution  might  exist  in  this  case. 


Value  of  Endowments 

Since  the  aggregate  demand  function  expresses  consumption  X<  =  Pt  •  Xt  as  a  function  of 
aggregate  income  1 t  =  Pt  •  It  and  prices,  we  reexpress  Theorem  1.5: 

st-Xt  =  (n’tst  -  It)s (KiPnt)-1  ■  H'xt ,  Pt .  i;  <  It  <  Pt  ■  /;* ,  (4.6.0) 

t 

where  Pt  is  total  size  of  the  population  in  period  t,  St  =  Pt  ■  St,  and  nt  =  8t~lft  are  discounted 
normalized  period  t  prices  ft]  aggregate  income  It  is  redefined  to  be  aggregate  income  measured 
in  terms  of  xt  instead  of  pt.  Before  we  place  Ft  (X* ) ,  the  inverse  of  this  demand  function,  into  the 
model,  we  must  first  relate  its  It  to  the  value  of  endowments.  This  is  important  because  the  model 
does  not  provide  any  detail  about  who  owns  the  endowments  and  therefore  no  detail  about  how 
rents,  wages,  dividends,  royalties,  interest  on  loans,  taxes,  government  doles,  etc.,  get  transferred  to 
the  final  consumers  for  consumption.  We  therefore  need  to  be  assured  that  the  model  nevertheless 
implicitly  provides  a  mechanism  whereby  the  total  value  of  endowments  used  for  consumption  are 
in  fact  transferred. 

Theorem  4.1.  The  value  of  endowments  It  used  to  produce  consumption  Xt  in  period  t  is  exactly 
equal  to  n^Xt,  the  attained  level  of  aggregate  income  used  for  consumption  in  period  t. 

Proof.  In  terms  of  the  prices  of  the  model,  the  value  of  endowments  available  to  period  t  is 
o't(Dt- iYt_i  -f  kt).  If  we  subtract  off  the  value  used  for  fixed  consumption  n[ ft ,  less  a{+1D(Yt 
the  value  passed  down  to  period  t  +  1,  the  net  by  definition  is  It: 

It  =  +  kt)  -  ft  —  cr't+iDtYt  (4.6.1) 

=  o'tBtYt  -  ffJ+iAYt  -  n'tft  (4.6.2) 

=  n't  ( At Y t  -  ft)  (4.6.3) 

where  (4.6.2)  follows  from  (4.1.0)  and  (4.5.0),  and  (4.6.3)  follows  from  (4.3.0)  and  (4.5.0).  Therefore 
by  (4.2.0)  and  (4.5.0): 

It  =  *t'Xt  .  |  (4.6.4) 

Thus,  independent  of  the  choice  of  the  demand  function,  the  primal  conditions  (4.1.0),  (4.2.0), 
the  profitability  conditions  (4.3.0)  and  the  complementarity  conditions  (4.5.0)  imply  that  It  =  x{Xt 
holds. 

Equivalent  Concave  Program 

Having  given  the  conditions  that  an  expected  per  capita  utility  function  and  inverse  of  the 
expected  per  capita  demand  function  exist,  we  now  assume  these  exist  and  set  aggregate  Xt  = 
PtXt,  St  =  Pt§,  Ut  =  Ptftt,  1 1  =  Pt  -  It  where  Pt  is  the  size  of  the  population  in  period  t.  Let 


p  =  G  1  (v)  where  v  =  §t  -  Xt 


(4.7.0) 


It  follows  from  (2.15.3),  after  rescaling  -l7t (.£*)  =  Z  for  population  size, 

U*(Xt)  =  -Pt  ■  Siip'&p)-'!*  (4.7.1) 

where  Qi  denotes  geometric  mean  of  (p'ff’p)-1/3  for  Hx ,  i  =  (1,. . .  ,n),  in  the  population  urn. 

Theorem  4.2.  The  equilibrium  problem  is  equivalent  to  solving  the  concave  program 

T 

maxU(X)  =  52^-1Ut(Xt)  (4.8.0) 

t=i 

subject  to  (Xt,  Yt)  >  0  and  the  primal  Bow  conditions  t  =  1, . . .  ,T: 

BtYt  <  Dt-iYt-i  +  kt  (4.8.1) 

—A* Yt  +  X<  <  ft  (4.8.2) 

providing  (i)  the  primal  problem  is  feasible,  (ii)  the  geometric  mean  of  (pJH*pt)1 12  fort  =  1,2,. ..,n 
is  used  as  a  numeraire  for  normalizing  prices  pt ,  and  (Hi)  the  aggregate  income  I(  =  xJXt  associated 
with  the  optimal  solution  satisBes  Pt  •  /*  <  I*  <  Pt  •  /**,  and  (iv)  Xt  >  0.  Under  these  conditions 
an  equilibrium  solution  exists. 

Proof.  The  condition  that  I(  is  bounded  between  certain  lower  and  upper  limits  is  our  way  of 
saying  that  the  demand  function  of  all  individuals  j  is  of  the  form  S3  —  X3  —  (jtSj  —  J)(x'iPjr)-1  • 
H3x,  which  would  not  be  the  case  if  their  budgets  I  were  extremely  low  or  extremely  high.  The 
equilibrium  problem  obviously  has  no  solution  if  the  primal  problem  is  infeasible.  Here  in  Part 
IV,  we  assume  conditions  of  Theorem  3.3  hold  so  that  a  utility  function  for  each  period  exists.  By 
definition  of  a  utility  function,  Ut(Xt)  are  homogeneous  functions  of  degree  1  in  v  =  St  -Xt  which 
are  strictly  concave  in  v  except  along  rays  with  v  =  0  as  origin.  The  concave  program  under  these 
conditions  has  a  finite  optimal  solution.  The  Kuhn- Tucker  conditions  for  optimality  turn  out  to  be 
the  same  as  the  dual  constraints  of  the  equilibrium  problem,  (4.3.0),  (4.4.0),  and  (4.5.0). 

IfXj  >  0,  then  in  the  equilibrium  model,  by  (2.7.1)  and  (2.15.3): 


xt  =  pt(xt)  =  au  ,/ax, 

(4.9.0) 

=  dZ/dv  =  Z  •  p 

(4.9.1) 

=p/MP'n'pyJ* 

(4.9.2) 

=  Pt/ Sitft&pt)117 . 

(4.9.3) 

The  denominator  SiiPtH'Pt)1^  may  be  viewed  as  a  numeraire  for  normalizing  period  t  prices  pt. 

Therefore  if  a  feasible  solution  to  the  primal  exists,  an  optimal  feasible  solution  exists  that 
satisfies  the  Kuhn-Tucker  optimality  conditions,  which  is  the  same  as  saying  an  equilibrium  solution 
exists.  The  strict  concavity  implies  that  the  values  of  X<  and  xt  are  unique;  those  of  Yt  and  ot 
need  not  be  unique.  | 

Comment.  The  above  way  of  normalizing  p  is  different  from  the  conventional  one  p/e'p  where 
g1  =  (1/m-"  1/m)  and  m  is  the  number  of  components  in  the  vector  p.  Later  on  we  will  present 
evidence  why  this  way  of  normalizing  pt  is  just  as  satisfactory  if  not  superior. 
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Theorem  4.3.  The  existence  of  a  utility  function  U(X)  depends  on  how  the  numeraire  for  period 
t  prices  used  by  the  investor/producers  in  the  rate-of-retum  formulas  is  defined. 


Proof.  It  is  sufficient  to  demonstrate  this  when  H*  =  R  for  all ».  In  this  case  G(xt)  =  {*'tR*t)  1  ■ 
Hxt  and  the  expected  aggregate  demand  function  (4.6.0)  simplifies  to: 


v  =  St -Xt  =  {xtSt  -  IJWRxt)-1  •  Rxt  . 


(4.10.0) 


Dropping  the  income  factor  and  rescaling  xt  as  p,  then  for  this  special  case: 

v  =  {p'Hp )_1  •  Rp  , 

p  =  (v’Mv)"1  •  Mv  ,  where  M  =  H*1  . 


(4.11.1) 

(4.11.2) 


By  (4.7.1)  and  (4.9.2),  -Ut  simplifies  in  this  case  to 


-Ut  =  Z  =  {p'Hp)-1/2  =  {v'Mv)1/2 
=  {($  -  Xt)'M{St  -  Xt)]1/2 


(4.11.3) 

(4.11.4) 


By  (4.9.0),  (4.9.1),  (4.11.3)  and  (4.11.2),  x  simplifies  to 


x  =  6*-1*  =  S'-1!  p  =  6*-1  Mv/ {v'Mv)1/2 


(4.11.5) 


It  is  easy  to  show 


{t'R*)1/2  =  1  .  (4.11.6) 

and  therefore  the  implied  numeraire  for  normalizing  prices  is  {ptHpt)1/2 ■ 

Because  It  =  x'tXt,  the  right-hand  side  of  (4.10.0)  is  a  homogeneous  function  in  x  of  degree 


0,  implying  that  the  inverse  demand  function  Ft(Xt),  that  expresses  xt  in  terms  of  Xt  =  Pt  •  Xt, 
while  not  unique,  can  be  determined  uniquely  except  for  a  scale  factor.  Note  in  Theorem  4.2 
we  choose  one  particular  way,  see  (4.9.3).  When  X<  >  0,  relation  (4.4.0)  becomes  an  equation 
Xt  =  5<_1Ft(Xt)  where  Ft(Xt)  =  ftt  is  defined  to  be  normalized  period  t  prices.  We  therefore 
must  choose  the  proportionality  factor  so  that  Ft  (Xt )  automatically  generates  normalized  prices  of 
period  t  before  discounting. 

The  numeraire  for  normalizing  prices,  however,  can  be  chosen  in  more  than  one  way.  If  investors 
calculate  their  rate  of  return  based  on  period  t  prices  pt  normalized  by  Xt  =  pt/{p'tRpt)1/2,  for  the 
special  case  H'  =  R  for  all  »,  then  for  this  per  capita  demand  function  simplification  (4.10.0),  the 
inverse  per  capita  demand  function  as  we  have  shown  can  be  stated  explicitly  and  is  proportional 
to  M(St  -  Xt)  or  M(S «  -  X*),  see  (4.11.2).  Under  this  definition  of  normalization,  Ft(X<)  satisfies 
for  X«  >  0,  M  =  R -1 


Ft{Xt)  = 


I(S«  -X*)'Af(S<  -  X*)1/3 


Af (St  -  Xt), 


(4.12.0) 
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which  is  the  same  as  (4.11.5)  before  discounting.  Note  in  this  case,  f'Hx  =  1  for  all  x  =  Ft(Xt) 
since  [Ft(Xt)'  H Ft  (X*]1/2  =  1  for  all  Xt .  This  way  of  normalizing  Ft  (Xt )  does  depend  on  the  scaling 
of  M  and  we  therefore  would  like  the  scaling  of  M  to  be  such  that  base-year  prices  p0  =  (1, . . . ,  l)' 
remain  the  same  before  and  after  normalization.  This  is  why  earlier,  we  required  H  to  satisfy 
1  =  p'0Hpo  =  Efc  E*  8kt,  see  (1.1.1). 

On  the  other  hand,  if  investors  calculate  their  rate  of  return  based  on  the  vector  of  prices  pt 
normalized  within  each  period  so  that  their  average  price  is  unity,  i.e.  normalized  by  xt  =  pt/e'pt 
where  e*  =  (1/m, . . . ,  1/m),  then  the  scale  factor  must  be  chosen  so  that  normalized  prices  satisfy 
e*x  —  1.  Ft  (Xt )  under  this  definition  of  normalization  satisfies  for  the  special  case  Hl  =  H  for  all 


Ft(Xt)  = 


M( St  -  Xt),  for  X*  >  0  ,  M  =  H 


-  rr-i 


(4.13.0) 


?M(St-X<)  v  ‘  .  '  ' 

Note  i? Ft{Xt)  =  1  for  all  Xt.  Moreover  under  this  definition,  Ft(Xt)  does  not  depend  on  the 
scaling  of  M.  Also  when  the  physical  units  for  categories  k  are  chosen  so  that  base  year  prices 
po  =  e  =  (1, . . . ,  1)',  they  remain  unchanged  when  we  normalize  in  this  way. 

For  the  H  estimated  by  Hu  Hui  from  empirical  data  [10]  and  Tabulated  in  Table  1,  the  normal¬ 
ized  prices  x<  —  Pt/ip'tHpt)1^  had  average  values  g'jrt  for  years  t  from  1961  to  1982  which  differed 
only  slightly  from  1,  see  the  last  column  of  Table  4,  so  that  the  investor  would  be  indifferent  as  to 
whether  the  denominator  of  (4.12.0)  or  (4.13.0)  were  used  for  numeraire. 

We  need  only  show  that  a  utility  function  does  not  exist  for  this  special  case  of  the  equilibrium 
problem  when  prices  are  normalized  by  Wt  =  pt/ept  and  F((Xt)  is  defined  by  (4.13.0).  Consider 
a  two  period  model  so  that  we  are  maximizing  the  utility  t/(X)  =  Ui(Xi)  +  6U2(X2)  subject  to 
the  primal  constraints  (4.1.0)  and  (4.2.0).  Further  suppose  Xi  has  only  two  components  so  that 
X*  =  (Xn,Xi2).  Let  Si  =  (Sn,S12),  Si  —  Xi  =  (Sn  —  Xn,Si2  —  Xi2).  Let  =  (xn,xi2) 
and  M  =  [m,-y]  be  a  2  x  2  symmetric  positive  definite  matrix.  Let  Si  -  Xi  =  V  =  (Vi,V2)  where 


Vi  —  Sn  —  Xu,  V2  —  S12  -  Xi2. 


(4.14.0) 


At  a  maximum  the  Kuhn-Tucker  conditions  dUi/dXi  —  xi  should  hold,  see  [12].  From  (4.13.0), 
x  =  MV jV  where  the  denominator  D  =  e'MV  —  |(mi;  -f  mi2)Vi  -(-  §(mJ2  +  m22)V2.  Therefore, 
the  following  should  hold: 


3Ui/dXn  —  xu  —  (mnVi  +  rni2V2)/P, 
dUi/5Xi2  =  xi2  =  (mi2Vi  +  m22V2)/P. 

In  order  for  a  utility  function  to  exist,  the  second  partial 

a2Ui/axn3Xi2  =  -d2v1/dxndv2 

computed  from  (4.15.1)  should  agree  with 

d2U1/dXl2dXu  =  -d2U1/dX12dVl 


(4.15.1) 

(4.15.2) 
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computed  from  (4.15.2)  for  a 11  choices  of  Xu, X12.  Setting  these  second  partials  equal  to  each 
other,  we  should  have 


n?i2  (mnVi  +  mi2V2)(mi2  +  m22)(l/2) 

D  P2 

-  -HU2  +  +  m22Vi)(mn  +  mi2)(l/2) 

"  D  P2 


(4.15.3) 


which  reduces  to  (mnm^j  -  m22)(Vri  -  V2)/P2  =  0.  But  this  does  not  hold  for  all  choices 
of  (Xn,Xi2)  because  M  being  positive  definite  implies  -  m22  >  0  and  (Vx  -  V2)  = 

(Su  -  Xu)  -  (S12  —  X12)  j*  0  for  general,  X,  a  contradiction.  | 

Our  proof  of  non-existence  of  an  objective  function  when  prices  are  normalized  in  the  conven¬ 
tional  way  is  only  valid  for  two  or  more  periods.  In  a  one  period  case,  the  factor  e'p  normalizing 
prices  in  the  rate  of  return  formula  can  be  factored  out  and  replaced  by  ( p'Hp )1^2.  The  new 
problem  is  then  equivalent  to  the  original  problem  but  in  its  new  form  a  utility  function  exists. 
However,  in  the  multiperiod  case,  it  is  not  possible  to  factor  out  &pt  and  replace  it  by  (pj  Hpt)1^2 
without  invalidating  the  rate-of-return  relation  (4.3.0). 

Comment . 

This  ends  the  theoretical  derivation  of  an  objective  function  for  the  economy.  We  conclude  that 
the  economy  will  grow  if  it  has  the  resources  and  technology  to  grow  and  if  it  pays  according  to 
the  aggregate  “utility”  function  to  trade  off  mo /erne  nt  of  the  consumption  vector  Xt  towards  the 
“satiation”  vector  St  of  earlier  periods  t  for  considerably  larger  movements  towards  the  satiation 
vectors  of  later  periods. 

Recalling  v  =  S<  —  Xt,  X<  =  Pt  •  Xt,  S t  =  Pt  St  and  It  =  Pth  where  Pt  is  population  size,  we 
have  by  Euler’s  Theorem  for  a  general  homogenous  utility  function  of  degree  1: 


U t  =  PfZ  =  Pf  {dZ/dv')v 
=  Pt  ■*'v  =  Pf(xlSt-It) 


(4.16.0) 

(4.16.1) 

(4.16.2) 


Thus  the  disutility  is  the  sum  of  the  discounted  additional  aggregate  income  (measured  in  normal¬ 
ized  period  t  prices)  needed  to  purchase  the  “satiation”  vector  over  various  periods  t  where  the  less 
additional  income  required  the  higher  the  “standard  of  living” . 

Ideally  normalized  prices  jft  should  have  the  property  that  a  unit  amount  of  income  should  enable 
each  individual  j  to  purchase  at  prices  goods  X}  whose  utility  to  j  is  unity.  In  a  certain  geometric 
mean  sense  this  is  true  if  jft  are  prices  p<  normalized  by  the  numeraire  GeometricMeani(ptHtpt)i /2. 
For  a  representative  set  of  individuals  1  with  utility  functions  Ul(S',H')  for  (S'yH')  in  the  urn 
and  income  P ,  let  —U\  =  [(S’  -  —  X’)]1/2  be  the  negative  utility  of  optimal  X1  to  «.  It 

is  easy  to  prove 

-Ul  =  (4.17.0) 
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where  J‘  =  f'tS%  —  /'  is  the  additional  income  »  requires  to  reach  satiation.  Letting  p,  denote 
geoemtric  mean,  it  follows  that  the  §*(—&{) / 5i(Jl)  is  unity  since  //’tt*)1/2  is  unity  when 
$i(ptHlpt J1/3  is  used  as  numeraire. 
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PART  V: 

SIMPLIFYING  AND  ESTIMATING  THE  DEMAND  AND  UTILITY  FUNCTIONS 

Simplifying  the  Form  of  the  Per  Capita  Demand  and  Utility  Functions. 

Our  immediate  goal  is  to  replace  the  key  factors  G(p)  =  £<[(p'#‘p)-1  •  H'p]  associated  with 
the  demand  function  and  -2  log  Z  =  £{  log(p'7/‘p)  associated  with  the  utility  function  by  simpler 
expressions  whose  parameters  are  easier  to  estimate  and  then  give  reasons  why  these  approximations 
may  be  very  good.  We  will  use  the  symbol  =  to  denote  approximately  equal. 

Note  that  the  demand  function  (1.6.0)  of  individual  i  is  invariant  to  the  scaling  of  Hl .  Positive 
definite  matrices  H'  —  (#**]  satisfy  e'H'e  =  ^kt  >  ®  where  e'  =  (1, . . . ,  1).  Therefore,  see 

(1.1.1),  we  have  assumed  without  loss  of  generality  that  H *  has  been  rescaled  so  that 


e'Hie  =  J2Y,Hkt  =  l  forall,‘>  e'  =  (l,...,l). 


(5.1.0) 


The  units  for  measuring  the  consumption  of  goods  in  “physical”  terms  are  usually  defined  so  that 
their  base  year  prices  per  unit  is  po  =  e  =  (1, . . . ,  1)'.  Therefore  p'0H'po  =  1  for  all  »  by  (5.1.0), 
implying  p'0Hpo  =  1  where  we  define  H  —  £tH‘ .  It  is  convenient  for  the  simplification  that  we  use 
as  numeraire  (p'Hp)1/7  implying  p’Hp  =  1  for  all  “normalized”  vectors  p.  It  follows,  in  particular, 
that  Po  =  Po/ip'oHpo)1^  =  Po  =  «  and  p^H'po  —  1  for  all  »  for  the  base  year. 

Theorem  5.1.  The  Brat-order  approximation  of  log  Z  as  a  function  of  p  —  G~1{St  -  Xt)  yields 
the  approximation: 

Z  =  9i  (pfTpy1'*  =  (P'HP)-1/2  ,  (5-2.0) 

where  positive  definite  Hx  are  rescaled  so  that  e'H'e  =  1  for  all  i  and  H  =  £iHx  and  Q,  denotes 
the  geometric  mean. 

Proof.  We  normalize  p  by  p  =  p/(p' Hp)1/2.  Subtracting  log (pHp)  from  both  sides  of  — 21og(Z), 
and  letting  p  =  e  +  A,  we  obtain 

-21og(Z)  -  log (p'ffp)  =  £i  log [(p'Hip)/{p'Hp)}  =  £t  log(p'ITp)  (5.3.1) 

=  £i  log((e  +  A)'  if  (e  -f  A)]  (5.3.2) 

=  £i  \og(e'Hxe  +  2A'tf*e  +  A'HXA)  (5.3.3) 

=  &log(l  + 2A'ITe  +  A'WA).  (5.3.4) 

As  a  first  order  approximation  of  a  natural  log,  set  log(l  +  e,)  =  e*. 

-21og(Z)  -  log {p'tip)  =  £,■[( ZA'We)  +  (A'tf’A)]  =  2A 'He  +  AHA  (5.4.1) 

=  (1  +  2A 'He  +  AHA)  -  1  =  (e  +  A )'H{e  +  A)  -  1  (5.4.2) 

=  ^ffp-  1  =  1- 1  =  0,  (5.4.3) 


where  p Hp  =  1  was  shown  in  comments  following  (5.1.0).  Therefore,  the  approximation  of  Z(p) 
is  given  by  (5.2.0).  | 

Little,  of  course,  is  known  about  the  utility  functions  of  individuals  or  how  they  vary  about  H. 
For  the  experiment  we  are  about  to  describe,  our  particular  choice  of  distribution  for  h,  =  p'  H'p,  is 
designed  to  show  that  the  approximations  of  Z(p)  and  hence  G(p)  =  -3  (log  Z)/dp  are  good  ones 
for  variable  p  in  the  neighborhood  of  some  fixed  po  =  (1, 1,  ...,)*  =  e  even  if  the  H'  for  individuals 
i  were  to  vary  a  lot.  In  Table  4,  prices  pt  for  the  years  t  =  1961  to  1982  relative  to  base  year  1972 
prices  p0  =  e,  are  tabulated.  For  p  in  the  neighborhood  of  p0  =  e,  the  values  of  h,  =  p  H'p  will 
obviously  be  distributed  around  the  mean  p1  Hp  =  1  with  a  very  low  standard  deviation,  because 
PqH'Po  —  1  for  all  *.  The  vector  of  prices  for  1982  differed  the  most  from  the  base  year  with  some 
prices  differing  by  more  than  30%.  Therefore  pt  for  t  =  1982  was  selected  as  an  extreme  case  of  pt 
for  our  illustrative  example. 

The  following  experiment  was  then  made.  The  symmetric  matrix,  denoted  H  in  Table  1,  was 
inverted  to  produce  an  M  and  a  thousand  cases  of  random  Ml  were  generated  by  independently 
varying  each  symmetric  pair  M'kt  about  by  (1  +  9)Mki  where  6  was  binomially  distributed 
with  mean  0,  and  standard  deviation  .1,  and  maximum  range  -.2  <  9  <.  2.  Three  of  the  thousand 
cases  were  dropped  because  M'  turned  out  not  to  be  positive  definite.  The  inverses  [M1]"1  =  H' 
were  then  computed  and  Hl  replaced  by  rescaled  Hx  so  that  e'H'e  =  1.  We  will  now  call  H,  the 
average  of  these  Hl  and  denote  by  p  =  pt / (p't H Pt)1  ^  where  pt  is  the  price  vector  for  year  1982. 
This  makes  the  mean  of  pH'p  =  1.  The  standard  deviation  of  pH'p  for  the  sample  of  997  cases 
turned  out  to  be  less  that  .05  or  5%.  For  our  analysis,  we  exaggerated  the  standard  deviation  of 
hi  —  p' H'p  to  be  10%  instead  of  5%  and  to  have  a  maximum  range  of  20%  about  their  mean  of  1. 
If  5%  instead  of  10%  were  used  in  Theorem  5.2  below,  the  error  e  would  have  been  negligible. 

Theorem  5.2.  Let  p  —  pt/(ptHpt).  IF  the  distribution  of  h,  =  ( p’H'p )  >  0  is  binomial  with 
mean  of  1,  standard  deviation  —  .1,  and  a  range  .8  <  p  H'p  <  1.2,  then  the  percent  error  in  the 
approximation  of  Z  by  (pJ/Tp*)-1/2  is  less  than  0.5%: 


-21og(i?)  =  log[(l  +  e)(p' Hp)}  where  |c|  <  0.5%  . 


(5.5.0) 


Proof.  We  seek  an  error  bound  for  e.  From  (5.5.0)  and  (5.2.0), 

log[(l  +  c)(p'ti?pt)]  =  £,[log(P; iPpt)]  ,  (5.6.1) 

log[(l  +  e)]  =  ^[log(p/ H'p)\  ,  p  =  p/{p'Hpf12  ,  (5.6.2) 

Our  task  is  to  prove  |e|  <  .005.  We  assumed  that  h,  =  pH'p,  roughly  speaking,  has  a  truncated 
normal  distribution  with  mean  of  1,  a  standard  deviation  less  than  .1,  and  .8  <  h,  <  1.2,  namely 
we  assumed  it  to  be  the  binomial  distribution  with  values  (.8,  .9,  1.0,  1.1,  1.2)  and  corresponding 
probabilities  (1/16,4/16,6/16,4/16, 1/16).  For  this  experiment,  we  therefore  have: 

log(l  +  e)  =  Jq  1o8  8  +  Jq  1o?  9  +  yg  1°8  10  +  ^  log  1.1  +  ^  log  1.2  ,  (5.7.0) 


.  .v-  y.ov.v,  ’.v.vva>v" 


with  an  error  |e|  =  .00505.  We  assert  without  proof  that  any  continuous  unimodal  distribution 
with  same  mean,  standard  deviation,  and  truncation  would  yield  about  the  same  error  |e|  <  .005. 
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Because  \og{Z)  where  Z  is  the  disutility  function  is  so  closely  approximated  by  -5  log (p’tBpt) 
for  reasonable  variability  of  p'tHlpt,  we  apply  the  approximation  to 
-5(log  Z)/dp  to  obtain  our  approximation  for  the  key  factor  G(p)  of  the  demand  function: 

G(p)  =  Siip'H^)-1  •  Wp  =  -d(log  Z)/dp  (5.8.1) 

=  :|d[log(p,J?p)]/3p  =  (p'f?p)-1  •  Hp  for  all  p  near  po  .  (5.8.2) 

Hence  the  approximations: 

v  =  (p'Bp )-1  •  Hp  ,  p  =  («'M«)_1  •  Mv  ,  (5.8.3) 

where  M  =  H*1  and  v  =  St  —  Xt.  We  now  substitute  these  approximations  into  various  demand 
function  theorems  and  summarize  them  here: 

Demand  functions  of  individuals  with  income  I  (Theorem  1.1): 

S’  -  X’  =  (p'S>  -  I)[p'H’p)~x  •  H’p  (5.9.0) 

where  /*  <  /  <  /**  =  pS’ . 


The  approximation  for  Theorem  1.3,  the  expected  demand  function  of  individuals  whose  income 

is  /: 

S-X1  =  {p'S  -  fyiip'&p)-1  •  iTp  =  ( p'S  -  I)(p'Hp)~l  •  Hp  (5.10.0) 

where  max  /*  =  /*</<  /**  =  min  p'S1,  Hl  rescaled  so  that  e'H'e  =  1  and  H  =  St  H' . 

The  approximation  for  Theorem  1.5,  the  expected  per  capita  demand  function: 

S-X  =  {p'S  -  ^^(p'fl-’p)-1  •  IPp  =  (p'S  -  I)[p'Bp r1  •  Hp  (5.11.0) 

where  /*</</**. 

The  approximation  for  Theorem  2.9,  the  utility  function  Ut  =  Ut{Xt)  for  period  t  which 
implies  the  per  capita  demand  function  (5.11.0)  is: 


Ut  =  -9i(p'Hip)~1/2  (5.12.1) 

=  -(; p'Hp)~ 1/2  =  -(v'Mv)1/2  (5.12.2) 

^-[(S-^y^S-X,)]1/2,  (5.12.3) 

where  M  =  B~l  and  Qi  denotes  the  geometric  mean  for  H'  in  the  population  urn.  We  now  apply 
the  approximation  (5.12.3)  to  Theorem  3.2,  after  rescaling  for  population  size. 
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Theorem  5.3.  The  equilibrium  problem  is  approximately  equivalent  to  solving  the  concave  pro¬ 


gram: 


max  U(X)  =  -  -  Xt)'Aft(St  -  X,)]1'3, 


(5.13.0) 


subject  to  (Xt,  Yt)  >  0  and  the  primal  Bow  conditions  t  = 

Bt Yt  <  Dt- iYt_i  +  kt  (5.13.1) 

-AtYt  +Xt  <  ft  (5.13.2) 

where  St,X(,Yt  are  aggregate  quantities,  and  M*  =  ,  Bt  =  C, H\ ,  and  H\  have  been  rescaled 

so  that  J2k  52 1  =  1. 

EMPIRICAL  EVIDENCE  OF  LINEARITY  OF  THE  DEMAND  FUNCTION  AT 
FIXED  PRICES 

M.  Avriel  in  his  studies  [2]  convoluted  the  income  distribution  with  average  observed  personal 
consumption  data  for  over  fifty  commodities  as  a  function  of  average  observed  income  per  person 
at  various  household  income  levels.  More  recently,  one  of  the  authors,  McAllister,  repeated  the 
same  experiment  with  more  recent  survey  data  for  certain  key  commodities.  For  their  studies  a 
distribution  of  income  for  future  periods  was  assumed  to  be  a  certain  known  function  of  attained 
per  capita  income,  namely  (a)  that  Cfc(I),  the  average  of  consumption  of  individuals  at  income 
I  at  fixed  prices  p  of  commodity  k  will  not  change  in  the  future;  (b)  that  per  capita  income  / 
may  increase  with  time;  and  (c)  that  the  distribution  of  income  about  I  will  retain  its  same  shape 
when  rescaled  proportional  to  /.  The  base  year  distribution  was  based  on  survey  data.  They  then 
determined  per  capita  consumption  as  a  function  of  per  capita  income  by  the  convolution  formulas 
(1.10.1)  and  (1.10.2).  The  resulting  per  capita  demand  functions  turned  out  to  be  remarkably 
linear  at  fixed  prices  over  a  wide  range  of  per  capita  income.  See  references  [9,  2,  21]. 

Tables  2A  and  2B  tabulate  for  survey  years  1972-73  and  1980-81  the  average  consumption  of 
U.S.  per  person  in  various  household  income  classes  of  Food,  Clothing,  Housing,  Housing  Oper¬ 
ations,  Transportation,  Recreation,  Personal  Care,  and  All  Other.  Each  of  these  8  categories  of 
consumption  k  was  plotted  against  the  average  income  used  for  consumption  per  person  in  each 
household  income  class.  See  Figures  1  to  8  immediately  following  Tables  2A  and  2B.  Before  plot¬ 
ting,  however,  the  average  consumption  per  person  in  dollars  in  each  category  was  divided  by  its 
survey  year  price  thereby  converting  the  units  for  measuring  the  amount  of  consumption  to  “phys¬ 
ical”  units.  In  order  to  exhibit  the  plots  for  the  two  survey  years  for  each  category  and  make  them 
on  the  same  graph  comparable,  the  average  income  per  person  in  each  household  income  class  for 
the  1980-81  survey  in  1980$  was  deflated  to  1972$.  Prices  used  for  converting  the  consumption  in 
survey  year  dollars  to  “physical”  units  can  be  found  in  Table  4.  The  units  for  each  category  k  are 
chosen  so  that  price  in  1972-73  for  each  k  is  unity.  Prices  for  other  years  are  deflated  to  1972$. 

Since  the  vector  of  prices  in  each  survey  year  is  fixed,  the  regression  of  average  consumption 
of  category  k  per  person  in  a  household  income  class  versus  average  income  per  person  should  be 
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linear  for  “the  range  of  income  of  interest”  according  to  Theorems  1.2  and  1.3  or  its  simplified 
form  for  estimating  its  parameters,  (5.10.0).  The  reader  is  encouraged  to  ignore  the  trend  lines  on 
Figures  1  to  8  (which  will  be  explained  later)  and  to  study  the  sixteen  sets  of  plotted  points  and 
judge  for  himself  how  linear  they  are  for  these  very  broad  categories  of  consumption. 

Tables  3  and  4  tabulate  per  capita  consumption  and  prices  by  category  in  the  years  1961  - 
1982.  We  found  the  family  household  data,  Tables  2A  and  2B,  somewhat  inconsistent  with  per 
capita  data  of  Table  3  possibly  because  they  came  from  different  types  of  surveys.  It  seemed  best 
to  use  the  graphs  of  Tables  2A  and  2B  only  for  the  purpose  of  extrapolating  a  guess  of  S.  Using 
this  guess,  the  data  of  Tables  3  and  4  were  then  used  to  estimate  a  positive  definite  H . 

Estimating  S.  According  to  Theorems  1.2  or  1.3,  the  expected  demand  function  of  persons 
at  various  income  levels  is  linear  in  income  up  to  a  level  that  some  individuals  can  buy  their 
satiation  vector.  If  we  assume  for  the  moment  that  this  is  true  empirically  and  the  income  level 
just  sufficient  to  buy  the  expected  satiation  vector  pS  is  known  (and  not  much  different  for  different 
individuals  »),  then  the  satiation  value  S(k)  for  the  i-th  category  could  be  found  by  reading  off  the 
ordinate  value  when  the  abscissa  is  equal  to  pS.  Since  the  price  vectors  for  the  two  survey  years 
are  slightly  different  their  expected  satiation  income  pS  in  1972$  could  be  different,  but  it  is  not 
too  unreasonable  to  assume  for  purposes  of  roughly  estimating  S  that  the  two  expected  satiation 
incomes  are  equal. 

With  this  rational  in  mind,  straight  lines  were  fitted  to  the  data  for  survey  years  1972-73  and 
1980-81.  The  two  straight  lines  shown  on  Figures  1  to  8  for  each  category  k  are  “eyeball”  fits  to  the 
data  with  greater  weight  given  to  the  high  end.  It  was  assumed  that  at  some  very  high  income  level 
the  ordinates  of  the  two  fitted  lines  for  each  category  k  would  be  sufficiently  close  to  one  another 
that  making  the  two  ordinate  values  equal  would  distort  very  little  the  fits  to  the  observed  data 
at  the  much  lower  income  levels.  We  arbitrarily  pegged  this  high  per  capita  income  level  to  be 
$25,510,  an  income  level  about  equivalent  to  that  of  a  three  person  household  income  of  $200,000 
in  1987  dollars.  Therefore,  the  abscissa  where  the  two  lines  intersect  is  at  $25,510.  Their  common 
ordinate  value  is  the  estimate  of  S[k).  This  was  done  for  each  k  except  for  k  =  Recreation  whose 
trend  line  in  1972-73  was  ignored  because  it  had  a  radically  higher  slope  from  that  of  1980-81.  On 
Figure  8  only  one  trend  line  is  shown;  this  means  that  two  trend  lines  were  estimated  to  be  the 
same.  The  estimates  for  S(k)  are  tabulated  in  Table  1. 

Admittedly  this  is  a  pretty  crude  “eye  ball”  way  to  estimate  5,  nevertheless  the  resulting  linear 
fits  appear  to  be  reasonably  good  in  most  cuses.  The  real  test,  however,  is  not  how  good  is  the 
estimate  of  S  but  how  good  is  the  per  capita  demand  function  (5.11.0),  found  using  this  estimate 
of  S  and  price  and  per  capita  consumption  data  in  various  years  to  estimate  H. 

Estimating  the  positive  definite  matrix  H  =  M~l . 

The  values  estimated  for  the  elements  of  matrix  H  are  also  tabulated  in  Table  1 .  The  matrix 
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H  =  [Si])  was  estimated  by  a  technique  of  solving  a  semi-infinite  linear  program  recently  developed 
by  Hu  Hui  [10]  that  finds  a  least  square  fit  under  the  restrictions  that  H  be  symmetric,  positive 
definite  and  EE Hij  =  1.  The  least  square  fit  was  made  to  (5.11.0)  after  multiplying  the  equation 
by  p'Hp.  The  source  data  for  this  estimate  was  (a)  the  guess  of  satiation  levels  S  described  above, 
and  (b)  the  observed  JCt  ,  which  we  denote  by  Xt ,  and  observed  prices  pt  from  the  national  income 
and  product  accounts  of  the  U.S.  for  years  t  from  1961  to  1982.  Per  capita  income  was  assumed 
to  be  the  observed  values  It  =  p'tXt-  Observed  Xt,  which  are  indices  of  physical  quantities  of 
consumption,  are  tabulated  in  1972$  in  Table  3  for  each  of  the  eight  categories  for  each  year  t  from 
1961  to  1982.  The  corresponding  prices  pt  for  these  items  deflated  to  base  year  1972  prices  are 
tabulated  in  Table  4.  The  predicted  values  of  Xt,  using  the  per  capita  demand  function  (5.11.0) 
with  the  estimated  parameter  values  given  in  Table  1,  are  tabulated  in  Table  5  for  each  of  the 
22  years.  A  comparison  of  how  well  predicted  Xt  compares  with  observed  Xt  can  be  made  by 
comparing  Table  3  with  Table  5  or  by  comparing  the  solid  line  curves  with  the  broken  line  curves 
in  Figures  9  to  16  immediately  following  Table  5.  Note  the  excellence  of  the  fit.  The  average 
percent  error  of  fit  can  be  found  at  the  bottom  of  T&ble  5. 

The  average  price,  after  the  price  vector  pt  is  normalized  by  using  (pJiZpt)1/2  as  numeraire, 
are  tabulated  in  the  last  column  of  Table  4.  It  can  be  seen  that  these  differ  only  slightly  from  unity, 
which  would  have  been  their  average,  had  they  been  normalized  using  average  price  as  numeraire. 
Thus  an  investor  would  be  quite  indifferent  as  to  which  of  the  two  ways  is  chosen  to  normalize 
future  prices  before  calculating  a  rate  of  return. 
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ESTIMATED  DEMAND  FUNCTION:  S  -  Xt  =  &L-J*.  .  Hpt 

PtHpt 


Entries  in  columns  below  are  estimated  Hkt  X  1000  for  k  >  £ 


k 

Est. 

Food 

*  =  (1) 

Cloth 

(2) 

Housing 

(3) 

House 

Oper. 

(4) 

Transp. 

(5) 

Recrea-  Per. 

tion  Care  Other 

(6)  (7)  (8) 

(1)  Food 

2666 

30 

(2)  Clothing 

1970 

1 

35 

(3)  Housing 

4670 

7 

33 

77 

(4)  House  Oper. 

2050 

-  6 

-14 

20 

110 

(5)  Transp. 

3479 

10 

9 

8 

54 

54 

(6)  Recreation 

2253 

18 

-17 

-  2 

-  3 

20 

56 

(7)  Per.  Care 

526 

12 

20 

-  7 

-77 

-  1 

-  1  97 

(8)  Other 

7340 

16 

15 

54 

49 

32 

23  -18  145 

Entries  Ski  above  the  diagonal  =  Htk  below  the  diagonal 

Pt 

Xt 

It 

S  = 

H 

Xt 


=  observed  price  vector  (period  t)  from  time  series  (Table  3) 

=  observed  per  capita  consumption  vector  from  time  series  (Table  4) 

=  observed  per  capita  toted  consumption  income  =  p[Xt  (Total  Column,  Table  4) 

=  estimated  per  capita  satiation  vector 
=  estimated  positive  definite  symmetric  matrix* 

=  per  capita  consumption  vector  predicted  by  demand  function 


*  Estimated  by  a  method  developed  by  Hu  Hui  that  yields  least  square  fit  subject  to  H  being  positive 
definite  [10], 
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I  liu,  biuimry  wus  prepared  by  Ujfotliy  Sltl t  leid;  Figures  1  lo  H  ilkit  follow  display  the  above  ^rajihicully.  Tliesje  were 
jncpaml  ty  Isabel  Pereira. 


Figure  1 


FOOD,  TOBACCO,  ALCOHOL  EXPENDITURES 

AS  A  FUNCTION  OF  PERSONAL  INCOME  FOR  CONSUMPTION 


o  =  1972-73  SURVEY 
•  =  1980-81  SURVEY  (DEFLAT0R=2.0375) 


isoo  aooe  ism  aooe  mm  som  «sm  som  ssoe  eooo 
TOTAL  INCOME  FOR  CONSUMPTION  (1978  DOLLARS) 
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Figure  2 


CLOTHING  EXPENDITURES 

AS  A  FUNCTION  OF  PERSONAL  INCOME  FOR  CONSUMPTION 


o  =  1972-73  SURVEY  / / 

•  =  1980-81  SURVEY  (DEFLAT0R=1.385T/  / 


isoo  aooe  asoo  aooo  asoo  «ooo  ««•  sooo  saoe  «ooo 
TOTAL  INCOME  FOR  CONSUMPTION  (1972  DOLLARS) 
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HOUSEHOLD  EXPENDITURES 

AS  A  FUNCTION  OF  PERSONAL  INCOME  FOR  CONSUMPTION 


tSOO  2000  2300  3000  3300  4000  4000  9000  3300  3000 

TOTAL  INCOME  FOR  CONSUMPTION  (1978  DOLLARS) 


Figure  6 


RECREATION  EXPENDITURES 

AS  A  FUNCTION  OF  PERSONAL  INCOME  FOR  CONSUMPTION 


°  =  1978-73  SURVEY 
•  =  1980-81  SURVEY  (DEFLATORsl.51935) 
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1500  3000  2300  3000  3500  4000  4300  5000  3500  5000 

TOTAL  INCOME  FOR  CONSUMPTION  (1972  DOLLARS) 


Figure- 7 

PERSONAL  CARE  EXPENDITURES 

AS  A  FUNCTION  OF  PERSONAL  INCOME  FOR  CONSUMPTION 
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Figure  8 

OTHER  EXPENDITURES 

AS  A  FUNCTION  OF  PERSONAL  INCOME  FOR  CONSUMPTION 


o  =  1972-73  SURVEY 
•  =  1980-81  SURVEY  (DEFLA.T0Rs2.08) 
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SOLRCZ:  National  inosm  and  Product  Accounts,  special  atypl— nt  to  the  Survey  of  Current 
Business  Statistic  Tahlaa  published  Sepcober  1961  for  yaara  1929-76,  1976-79  and 
Ravlsad  Eatlxaatea  published  July  1983/Vol.  63,  No*  7.  Table  2.4  Personal 
Consmpcion  tlxpendi cures  by  Type  of  E^aendlture. 


Table  4 

p  OBSERVED  FUCKS  ST  THAIS  RELATIVE  TO  BASK  TUI  1972 
C  BT  UFA  CAIBCOUKS 


YEAR 

FOOD 

CLOTH 

HOUSING 

HOUSE 

0PER. 

TRANSP. 

RECREA¬ 

TION 

PER. 

CARE 

OTHER 

eii 

1961 

.97 

1.01 

1.07 

1.05 

1.08 

1.07 

1.03 

.86 

1.03 

1962 

.96 

1.00 

1.06 

1.04 

1.08 

1.07 

1.02 

.88 

1.02 

1963 

.96 

1.00 

1.06 

1.06 

1.04 

1.07 

1.03 

.91 

1.02 

1964 

.97 

.99 

1.05 

1.04 

1.06 

1.07 

1.02 

.89 

1.02 

1963 

.97 

.98 

1.04 

1.03 

1.06 

1.07 

1.03 

.92 

1.01 

1966 

.99 

.98 

1.03 

1.03 

1.05 

1.05 

1.03 

.92 

1.01 

1967 

.98 

1.00 

1.02 

1.02 

1.04 

1.04 

1.03 

.93 

1.02 

1968 

.98 

1.02 

1.01 

1.03 

1.03 

1.04 

1.04 

.96 

1.01 

1969 

.99 

1.03 

1.00 

1.02 

1.02 

1.03 

1.02 

.97 

1.01 

1970 

1.00 

1.02 

.99 

1.02 

1.01 

1.01 

1.02 

.98 

1.01 

1971 

.99 

1.01 

1.00 

1.01 

1.02 

1.00 

1.01 

.97 

1.00 

*1972 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1973 

1.05 

.98 

1.00 

.97 

.97 

.97 

.99 

1.01 

1.00 

1974 

1.08 

.95 

.98 

.98 

.99 

.93 

.99 

.99 

1.00 

1975 

1.07 

.91 

.98 

1.01 

.99 

.91 

1.03 

1.00 

1.00 

1976 

1.05 

.90 

.99 

1.03 

1.01 

.90 

1.05 

1.03 

.99 

1977 

1.04 

.88 

1.00 

1.00 

1.03 

.88 

1.06 

1.04 

.99 

1978 

1.06 

.84 

1.00 

1.01 

1.02 

.85 

1.07 

1.06 

.99 

1979 

1.07 

.80 

1.01 

.99 

1.06 

.81 

1.07 

1.07 

.98 

1980 

1.05 

.75 

1.02 

.97 

1.12 

.79 

1.06 

1.07 

.97 

1981 

1.05 

.71 

1.03 

.97 

1.13 

.77 

1.06 

1.07 

.97 

1982 

1.04 

.69 

1.05 

.99 

1.09 

.77 

1.10 

1.08 

.97 

*  base  year 

SOURCE:  National  Income  and  Produce  Accounts,  special  supplement  to  the  Survey 
of  Currane  Buainaaa  Statistic  Tables  published  September  1981  for 
years  1929*76,  1976*79  and  Revised  Estimates  published  July  L983/Vol. 
63,  No*  7*  Table  2.4  Personal  Consumption  Expenditures  by  Type  of 
Expenditure*  Price  index  Is  derived  by  comparing  expenditures  in 
current  dollars  in  Table  2.5  with  those  in  Table  2.4* 

_  *  1/2 

eit£  -  average  value  of  components  of  i*  •  Pc/(p’Hpc> 
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Table  5 


f  PREDICTED  OOKSUMPTIOa  PQ  CAPITA  R  YEAIS  (1972$) 
e  IT  MIPA  CATBGOKIXS  IMP  OH  ODttMD  FODCTIOR 


YEAR 

FOOD 

CLOTH 

HOUSING 

HOUSE 

OPER. 

TRANSP. 

RECREATION 

PER. 

CARE 

1961 

703 

177 

431 

228 

322 

189 

68 

1962 

713 

194 

454 

219 

327 

185 

89 

1963 

723 

200 

466 

223 

351 

196 

85 

1964 

718 

226 

498 

227 

362 

181 

97 

1963 

728 

237 

534 

255 

382 

194 

77 

1966 

731 

242 

558 

267 

408 

219 

69 

1967 

744 

228 

562 

295 

425 

245 

47 

1968 

763 

236 

587 

293 

456 

271 

52 

1969 

772 

242 

601 

288 

475 

287 

75 

1970 

775 

249 

615 

281 

492 

304 

78 

1971 

785 

246 

616 

286 

493 

315 

86 

1972 

798 

271 

645 

283 

534 

325 

110 

1973 

783 

273 

663 

337 

578 

342 

77 

1974 

767 

282 

682 

312 

556 

346 

81 

1973 

778 

299 

703 

283 

573 

363 

71 

1976 

812 

314 

735 

260 

595 

397 

96 

1977 

835 

317 

763 

291 

610 

423 

80 

1978 

839 

334 

799 

312 

645 

442 

64 

1979 

845 

335 

821 

334 

636 

459 

54 

1980 

852 

344 

838 

329 

592 

444 

72 

1981 

859 

359 

850 

311 

599 

451 

90 

1982 

864 

357 

837 

313 

643 

469 

57 

S 

2666 

1970 

4670 

2050 

3479 

2253 

Avg.  Z 
| error | 

.9 

3.4 

1.7 

3.4 

2.1 

3.8 

Z  | error)  •  100  are  tabulated  la  Table  3 

Xt  computed  using  per  capita  demand  function  given  In  Table  l 

Graphe  that  follow  thla  Table  conparlng  St  with  St  are 
by  Isabel  Pereira. 


OTHER 


*34 
443 
410 
477 
485 
523 
536 
543 
556 
558 
585 
598 
633 
64  L 
640 
659 
706 
723 
7  37 
754 
766 
760 
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Figure  9 


FOOD.  TOBACCO.  ALCOHOL  EXPENDITURES  PER  CAPITA 

FOR  YEARS  1961-1982 
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Figure  10 


CLOTHING  EXPENDITURES  PER  CAPITA 

rOR  YEARS  1961-1982 
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Figure  11 


HOUSING  EXPENDITURES  PER  CAPITA 

FOR  YEARS  1961-1902 
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Figure  12 
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Figure  13 


TRANSPORTATION  EXPENDITURES  PER  CAPITA 

FOR  YEARS  1961-1982 
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PERSONAL  CARE  EXPENDITURES  PER  CAPITA 

FOR  YEARS  1961-1982 
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Figure  10 


OTHER  EXPENDITURES  PER  CAPITA 

FOR  YEARS  1961-1902 
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ABSTRACT 

Given  a  general  dynamic  equilibrium  formulation  of  a  time  staged  model,  we  seek  conditions  on  the 
distribution  of  utility  functions  of  individuals  which  imply  the  model  is  equivalent  to  a  mathematical  program. 

Gorman  and  others  long  ago  have  observed  that  Engel  curves  of  average  consumption  as  a  function 
of  income  at  fixed  prices  are  remarkably  linear  over  a  broad  range  of  income  of  interest  which  tapers  off 
at  both  ends  of  this  range.  We  reproduce  this  phenomenon  by  assuming  (a)  that  a  general  polynomial  of 
the  second  degree  has  enough  parameters  (coefficients)  to  globally  represent  the  utility  functions  of 
individual  consumers,  and  (b)  the  distribution  of  utility  functions  that  individuals  have  is  independent  of  the 
income  they  happen  to  have.  We  achieve  the  latter  by  assigning  values  to  the  parameters  of  the  utility 
functions  by  a  random  drawing  with  replacement  from  a  'population  urn*  containing  a  representative  sets 
of  the  parameters.  We  then  derive  the  functional  form  of  the  per  capita  demand  function  and  necessary 
and  sufficient  conditions  for  its  integrability. 

Finally,  we  show  in  the  context  of  the  time  staged  model,  that  when  the  population  is  not  too 
polarized  as  to  its  tastes  at  fixed  income  levels,  a  concave  objective  function  always  exists,  which 
maximized  subject  to  the  physical  flow  constraints,  implies  the  equilibrium  conditions. 
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